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We investigate the dynamical and chemical evolution of the Large Magellanic 
Cloud (LMC) interacting with the Galaxy and the Small Magellanic Cloud (SMC) 
based on a series of self-consistent chemodynamical simulations. Our numerical models 
are aimed at explaining the entire properties of the LMC, i.e., the observed structure 
^ I and kinematics of its stellar halo and disk components as well as the populations 

. of the field stars and star clusters. The main results of the present simulations are 

summarized as follows. 

(1) Tidal interaction between the Clouds and the Galaxy during the last 9 Gyr 
transforms the initially thin, non-barred LMC disk into the three different components: 
the central bar, thick disk, and kinematically hot stellar halo. The central bar is 

, composed both of old field stars and newly formed ones with each fraction being equal 

' in its innermost part. The final thick disk has the central velocity dispersion of ~ 30 

' km s~^ and shows rotationally supported kinematics with Vin/ao ~ 2.3. 

(2) The stellar halo is formed during the interaction, consisting mainly of old stars 
Q ■ originating from the outer part of the initially thin LMC disk. The outer halo shows 
|H \ velocity dispersion of ^ 40 km s^^ at the distance of 7.5 kpc from the LMC center and 
C/3 . has somewhat inhomogeneous distribution of stars. The stellar halo contains relatively 

' young, metal- rich stars with the mass fraction of 2 %. 

(3) Repetitive interaction between the Clouds and the Galaxy enhances moder- 
. ately the star formation rate to ~ 0.4 Mq yr^^ in the LMC disk. Most of the new 

^ ' stars (~ 90 %) are formed within the central 3 kpc of the disk, in particular, within 

the central bar for the last 9 Gyr. Consequently, the half mass radius is different by a 
factor of 2.3 between old field stars and newly formed ones. 

(4) Efficient globular cluster formation does not occur until the LMC starts inter- 
acting violently and closely with the SMC 3 Gyrs ago). The newly formed globular 
cluster system has a disky distribution with rotational kinematics and its mean metal- 
licity is ~ 1.2 higher than that of new field stars because of the pre-enrichment by the 
formation of field stars prior to cluster formation. 

(5) The LMC evolution depends on its initial mass and orbit with respect to the 
Galaxy and the SMC. In particular, the epoch of the bar and thick disk formation 
and the mass fraction of the stellar halo depend on the initial mass of the LMC. 

Based on these results, we discuss the entire formation history of the LMC, the 
possible fossil records of past interaction between the Clouds and the Galaxy, and the 
star formation history of the SMC for the last several Gyr. 

Key words: Magellanic Clouds - galaxies:structure - galaxies:kinematics and dy- 
namics - galaxies :halos - galaxies:star clusters 



1 INTRODUCTION 

Galaxy interaction is generally considered to play a ma- 
jor role not only in controlling formation histories of field 
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stars and star clusters (e.g., Kennicutt 1998; Ashman & 
Zepf 1992) but also in transforming galactic morphologies 
(e.g., Noguchi 1987). The Large and Small Magellanic Cloud 
(LMC and SMC), which are believed to be interacting with 
each other, have long boon served as an ideal laboratory to 
study the detail of the tidal effects on structural, kinemati- 
cal, and chemical properties of galaxies based on the compar- 
ison between observations and numerical simulations (e.g., 
Murai & Fujimoto 1980, hereafter MF; Gardiner, Sawa, 
& Fujimoto 1994, GSF; Gardiner & Noguchi 1996, GN; 
Yoshizawa & Noguchi 2003, YN). Most of previous theo- 
retical and numerical papers on tidal interaction between 
the Clouds and the Galaxy however have discussed the ori- 
gin of Magellanic stream and the evolution of the SMC (e.g., 
Lin & Lynden-BcU 1977, 1982; Mathewson et al. 1987; MF; 
GSF; GN; YN) rather than the formation and the evolution 
of the LMC itself. Given the fact that recent observations 
have raised and confirmed several important questions re- 
lated to the star formation history, the physical properties 
of star clusters, and the dynamical properties of the LMC, 
it is doubtlessly worthwhile to discuss a theoretical model 
which provides an integrated and systematic understanding 
of the formation and the evolution of the LMC (Westerlund 
1997). 

One of the long-standing and remarkable problems re- 
lated to the star formation history of field stars in the LMC 
is as to whether there were epochs of dramatic increase in 
star formation rate of the LMC disk a few or several Gyr 
ago (e.g.. Butcher 1977; Stryker 1983; Bertelfi et al. 1992; 
Olszewski et al. 1996; Gallagher et al. 1996; Vallenari et al. 
1996; Holzman et al. 1999). Several authors have recently 
investigated star formation histories of the different regions 
of the LMC disk and bar based on the color magnitude di- 
agrams of the field stars derived by Hubble Space Telescope 
(Ardeberg et al. 1997; Elston et al. 1997; Holzman et al. 
1999; Olsen 1999; Smeckcr-Hane et al. 2002), and revealed 
that the epoch of enhanced star formation and the degree of 
the enhancement arc different between different regions. For 
example, Smecker-Hane et al. (2002) found that (1) stellar 
populations within the LMC bar were formed in the episodes 
of star formation about 4 — 6 and 1-2 Gyr ago, and (2) these 
burst populations can account for ^ 25 % and ~ 15 % of the 
LMC's stellar mass, respectively. The origin of the observed 
spatially different stax formation histories in the LMC disk 
is one of key questions related to the LMC evolution (e.g., 
van den Bergh 2000a). 

Several important physical properties of the globular 
clusters and populous young blue clusters in the LMC are 
in a stark contrast to those in the Galaxy (e.g., van den 
Bergh 2000a). These include the more flattened shapes of 
the LMC clusters (e.g., Gcislcr & Hodge 1980; van den Bergh 
& Morbey 1984), the disky distribution of its globular clus- 
ter system (e.g., Schommer et al. 1992), possible rotational 
kinematics of old clusters (e.g.. Freeman et al. 1983), a larger 
fraction of apparently binary clusters or physical cluster 
pairs in the LMC (Bhatia & Hatzidimitriou 1988; Bhatia et 
al. 1991; Dicbafl & Grebel 1998), a possible "age/metallicity 
gap" (e.g.. Da Costa 1991; Olszewski et al. 1991; Geisler et 
al. 1997; Sarajedini 1998; Rich ct al. 2001), and larger sizes 
at a given galactocentric distance (van den Bergh 2000b). It 
is unclear whether these differences are understood in terms 
of the LMC being much more strongly influenced dynam- 



ically by other nearby galaxies (i.e., the Galaxy and the 
SMC) compared with the Galaxy. 

Structural and kinematical properties of the LMC have 
been investigated by many authors concerning different stel- 
lar populations and gaseous components (Hartwick & Cow- 
ley 1988; Meatheringham et al. 1988; Irwin 1991; Luks & 
Rohlfs 1992; Kunkel et al. 1997; Graff et al. 2000; Olsen & 
Salyk 2002; Cioni & Habing 2003; Staveley-smith et al. 2003; 
Subramaniam 2004; See Westerlund 1997 for a review). For 
example, Caldwell & Coulson (1986) found that the east side 
of the LMC is closer to us than the west based on photo- 
metric observations of carbon stars. Wide-field photometric 
observations of the LMC disk revealed the exponential scale 
length of 1.5 kpc (Bothun & Thompson 1988) whereas an 
exponential disk model that fits to the distribution of field 
RR Lyrae stars has a scale length of 2.6 kpc (Kinman et al. 
1991). Based on observational data of carbon stars by the 
Deep Near-Infrared Southern Sky Survey (DENIS) and Two 
Micron All-Sky Survey (2MASS), van den Marcl et al. (2002) 
(hereafter vdMAHS) have recently shown that the LMC has 
a considerable vertical thickness with V/a of 2.9±0.9 and 
a mass of (8.7±4.3)x 10^ M© within 8.9 kpc. The origin 
of the observed differences in dynamical properties between 
various populations in the LMC remains unclear, though 
such differences may well provide some information about 
physical roles of galaxy interaction in the LMC evolution 
(Westerlund 1997; van den Bergh 2000a). 

Only some theoretical attempts have been made to un- 
derstand these observed properties of the LMC disk (Pagel 
et al. 1998; GSF; Weinberg 2000; Kumai et al. 1993), while 
a growing number of observational results have emerged. 
Weinberg (2000) first investigated the dynamical effects of 
the Galactic tidal field on structure and kinematics of the 
LMC disk. Bekki et al. (2004a) demonstrated that a single 
or binary star cluster is formed in a cloud-cloud collision 
triggered by tidal interaction between the LMC and SMC. 
However, these studies have the following two disadvantages 
in understanding the LMC evolution in a comprehensive 
manner. Firstly, these models are not fully self-consistent 
as they consider either dynamical evolution alone or empir- 
ical one-zone chemical evolution. Secondly, the model pa- 
rameters adopted in these works did not account for the 
revised knowledge on the structure and kinematics of the 
LMC as revealed by the latest observations such as DENIS 
and 2MASS (van der Marel & Cioni 2001; vdMAHS). Thus, 
the integrated and systematic understanding of structural, 
kinematical, and chemical properties of the LMC is therefore 
yet to be obtained. 

Thus, the purpose of this paper is to investigate these 
unresolved problems on the LMC formation and evolution 
based on chemodynamical simulations with the model pa- 
rameters (such as LMC's mass) consistent with the latest 
observations. In particular, we examine the following five is- 
sues: (1) how the star formation history of the LMC is infiu- 
cnccd by tidal interaction between the LMC, the SMC, and 
the Galaxy, (2) whether the stellar halo formation process of 
the LMC is similar to or different from that of other mem- 
bers of the Local Group of galaxies (e.g., M33, the Galaxy, 
NGC 3109, and NGC 6822), (3) how we can explain the ori- 
gin of structure and kinematics of the thin/thick stellar disks 
of the LMC, (4) how the LMC's chemical evolution process 
is associated with the dynamical evolution influenced both 
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Figure 1. Schematic view of tlic Galaxy and tlie Magellanic 
Clouds. Tlic present position of the LMC and that of the SMC are 
indicated by a large filled circle and by a small one, respectively, 
in the standard galactocentric coordinate system (X, y, Z). The 
orbital evolution of the LMC and that of the SMC for the last 1 
Gyr (— 1 < T < Gyr) are shown by a solid line and by a dotted 
one, respectively. 



by the Galaxy and the SMC, and (5) whether the origin of 
the unique nature of the globular clusters and blue popu- 
lous clusters is understood in terms of the tidal interaction 
between the Clouds and the Galaxy. Based on these investi- 
gations, we attempt to provide a entire history of the LMC 
for the last 10 Gyr. 

The layout of this paper is as follows. In §2, we sum- 
marize our numerical models used in the present study and 
describe the methods for analyzing structure and kinematics 
of the simulated LMC. In §3, we present numerical results on 
the time evolution of morphology, metallicity distribution, 
and dynamical properties of the LMC. In §4, we discuss the 
above five outstanding issues related to formation and evo- 
lution of the LMC. The conclusions of the preset study are 
given in §5. 



2 MODEL 

The present investigation is two- fold. First, we derive the 
most plausible and realistic orbits of the Clouds with re- 
spect to the Galaxy by using a backward integration scheme 
(MF; GSF; GN). Then we investigate chemodynamical evo- 
lution of the LMC on the derived orbit by using a fully 
self-consistent N-body model of galaxies (Bekki & Shioya 
1998, 1999; Bekki & Chiba 2000, 2001). Since our main pur- 
pose of this paper is to discuss the origin of structural and 
kinematical properties of the LMC, wc briefly describe the 
results of the orbital evolution derived from the backward 
integration scheme in this section. We describe the results 
of the N-body simulations in the next section §3. 



2.1 Derivation of the orbits of the Clouds 

2.1.1 The backward integration scheme 

We adopt the backward integration scheme originally de- 
vised by MF to derive the most reasonable and realistic three 
dimensional (3D) orbits of the LMC and the SMC around 
the Galaxy. In calculating the orbital evolution of this triple 
interacting system, wc adopt the most recently derived ob- 
servational parameters on locations and radial velocities of 
the Clouds with respect to the Galaxy and masses of the 
Clouds. In order to derive the orbits of the Clouds, we need 
to assume model parameters for the following quantities of 
the Galaxy and the Clouds: (1) the shape of the Galactic 
potential as a function of the distance r from the Galac- 
tic center, in particular, beyond 200 kpc where the Clouds 
reach at their apocenter passages, (2) gravitational poten- 
tial of the Clouds, (3) total masses (or mass profiles) of the 
Clouds, (4) dynamical friction between the Galactic dark 
halo and the LMC (SMC), and (5) dynamical friction be- 
tween the LMC and the SMC. The above (3) and (5) are 
more carefully considered in the preset study as described 
below. 

The gravitational potential of the Galaxy 3>g is assumed 
to have the logarithmic potential; 

$G(r) = -Fo'lnr, (1) 

where r and Vq are the distance from the Galactic center 

and the constant rotational velocity (= 220 km s^^), respec- 
tively. The LMC is assumed to have the Plummer potential; 

'I>L(rL) = -MLMc/(rL' + aif\ (2) 

where Mlmc, '"l, and ol are the total mass of the LMC, 
the distance from the LMC, and the effective radius, respec- 
tively. We adopt the same value of (=3 kpc) as previ- 
ous numerical studies adopted (e.g., GN). Recent dynamical 
study of the LMC by vdMAHS showed that the dynamical 
mass within 8.9 kpc of the LMC is (8.7 ± 4.3) x 10® M©, 
which is less than the half of the mass (2.0 X 10^° Mq) 
adopted in previous numerical studies (e.g., GN). Consider- 
ing this latest and more robust estimation of the dynamical 
mass of the LMC, we mainly investigate the orbital models 
with Mlmc = 10^° Mq. The SMC is assumed to have the 
Plummer potential; 

$s(rL) = -GMsMc/(rs' + al)°\ (3) 

where Msmc, ^s, and as are the total mass of the SMC, 
the distance from the SMC, and the effective radius, respec- 
tively. We adopt the same values of ol (=3 kpc) and Msmc 
(=3.0 X 10® Mo) as previous numerical studies adopted 
(e.g., GN). 

We consider the dynamical friction due to the presence 
of the Galactic dark matter halo both for the LMC-Galaxy 
interaction and for the SMC-Galaxy one and adopt the fol- 
lowing expression (Binney & Tremaine 1987); 

FfricG = -0.428 In Ag , (4) 

where r is the distance of the LMC (the SMC) from the cen- 
ter of the Galaxy. The mass M is either Mlmc or Msmc, de- 
pending on which Cloud's orbit (i.e., LMC or SMC) we cal- 
culate. We adopt the reasonable value of 3.0 for the Coulomb 
logarithm Ag (GSF; GN) both in the orbital calculation of 
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Table 1. Orbit models 

'Mq) Dynamical friction (C/l, Vl, VI/l) (kms-^) (f7s, Vs, VKs) (kms-^) 

yes (-5,-225,194) (40,-185,171) 

no (-5,-225,194) (40,-185,171) 

yes (-5,-225,194) (40,-185,171) 

no (-5,-225,194) (40,-185,171) 

yes (-56,-219,186) (60,-174,173) 

yes (41,-200,169) (60,-174,173) 

yes (-15,-225,194) (40,-185,171) 

yes (5,-225,194) (40,-185,171) 

yes (-2,-225,194) (40,-185,171) 

yes (-5,-225,194) (50,-185,171) 

(-5,-225,194) (40,-185,171) 



model 


Mt ( XlO 


A 


1.0 


B 


1.0 


c 


2.0 


D 


2.0 


E 


1.0 


F 


1.0 
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1.0 


H 


1.0 


I 


1.0 


J 


1.0 


K 


1.0 



the LMC and in that of the SMC. The above equation (4) is 
essentially the same as that shown in the equation (18) by 
MF. In addition to the above Ffric.G, the dynamical friction 
between the LMC and the SMC is also considered in the 
present study. We adopt the following expression; 



JfricLS = -0.428 In Als 



GMs^ 



(5) 



where tls and Als are the distance between the center of 
the LMC and that of the SMC and the Coulomb logarithm, 
respectively. Jfric.LS is assumed to act on the SMC, only 
when the SMC is within the LMC's tidal radius rt within 
which dark matter halo is gravitationally bound without 
being stripped from the Galactic tidal field. By using the 
theoretical model adopted in the equation (4) of GN, n is 
estimated as 1,3 kpc for the present model. vdMAHS obser- 
vationally estimated rt as 15.0 ± 4.5 kpc based on the newly 
derived total mass of the LMC. Therefore our choice of rt = 
13 kpc is regarded as a quite reasonable value. For compar- 
ison, we investigate the models with rt = kpc, in which 
dynamical friction between the Clouds is not included at all. 

By integrating equations of the motions of the Clouds 
toward the past from the present epoch, we investigate 
orbital evolution of the Clouds for given initial positions 
and velocities of the Clouds. We adopt the reasonable sets 
of orbital parameters that are consistent with observa- 
tions and thus were adopted in previous numerical stud- 
ies (GN). Figure 1 shows a schematic view of the Galaxy 
and the Magellanic Clouds and the orbital evolution of the 
Clouds. The current Galactic coordinate {b,l), where I and 
b are the Galactic longitude and latitude, respectively, is 
(-32.89, 280.46) for the LMC and (-44.30, 302.79) for the 
SMC, and accordingly the current positions {X, Y, Z) in 
tmits of kpc in the figure are (—1.0, —40.8, —26.8) for the 
LMC and (13.6,-34.3,-39.8) for the SMC. The current 
distance and the Galactocentric radial velocity of the LMC 
(SMC) is 80 (7) km s-^ 

We try to investigate the LMC evolution between the 
epoch immediately after its thin disk formation and the 
present-day. Oswalt et al. (1996) estimated the oldest stellar 
components of the Galactic thin disk as 9.5 ± 1 Gyr based 
on the faint end of luminosity function of white dwarfs in the 
Galaxy. There are no observational studies that are based on 
the luminosity function of the white dwarfs of the LMC disk 
and thus reveal how old the oldest stellar components of the 
LMC disk axe. Therefore, as a compromise, we adopt the 
assumption that the age of the oldest stellar populations of 




Time (Gyr) 

Figure 2. The orbital evolution of the Clouds for the last ~ 9 Gyr 
(-8.8 < T < Gyr) for the best orbital model A. The distance 
between the Cloud, that between the Galaxy and the LMC, and 
that between the Galaxy and the SMC are represented by a thick 
solid line, a thin solid one, and a dotted one, respectively. 



the LMC is similar to that of the oldest ones in the Galaxy: 
We thus investigate the LMC evolution for the last 9 Gyr. 



2.1.2 The representative orbits 

The current space velocities or (U, V, W) in units of km 
s^^ arc the most important parameters that determine 
the orbital evolution of the Clouds in the present models. 
They are represented by {Ul,,Vl,Wl,) for the LMC and 
by (C^s,^s,W^s) for the SMC. Eleven representative mod- 
els with different values of (?7l, VL, T4^l) and {Us,Vs,Ws) 
are discussed in the present study, and the Table 1 summa- 
rizes the model parameters for these: Model number (col- 
umn 1), total mass of the LMC represented by Mlmc in 
units of 10^°Mq (2), whether or not the dynamical friction 
between the Clouds is included (3), ((7l,Vl,Wl) (4), and 
{Us,Vs,Ws) (5). 

Among these 11 models labeled as A— K, the model 
A with (C7l,Fl, Wl) = (-5,-225,194) and {Us,Vs,Ws) = 
(40, —185, 171) is considered to be the best orbital model 
and thus referred to as the "best orbital model" through- 
out this paper. This is mainly because the Magellanic 
stream is self-consistently reproduced in previous models for 
these values. Different observations suggested different val- 
ues of (t/L,li,WL) and {Us,Vs,Ws) (e.g., Kroupa & Bas- 
tian 1997, KB; vdMAHS). For example, ({/l.Vl.Wl) and 
{Us,Vs,Ws) are (41±44, -200±31, 169±37) km s'^ and 
(60±172, -174±172, 173±128) km s'S respectively, for KB. 
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Without dynamical friction 




Time (Gyr) 

Figure 3. The same as Figure 2 but for the orbital model C 
(lower) and D (upper). It is clear from this figures that the Clouds 
can not be in a binary state more than ^ 6 Gyr if the dynamical 
friction of the Clouds is included in the orbital calculation. 
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Figure 4. The same as Figure 2 but for the orbital model B 
(upper left), I (upper right), J (lower left), and F (lower right). 
The values of the model parameters for these orbits are shown in 
the Table 1. 



([/L,Vi,W^L) invdMAHS is (-56 ±36, -219 ±23, 186±35). 
The best model accordingly is broadly consistent with the 
above latest proper motion data (KB). Taking this differ- 
ence into account, we investigate the orbital models with 
the initial velocity components adopted from the above ob- 
servations: Model E from vdMAHS, model F from KB. We 
also investigate the orbital models with the initial velocity 
components slightly different from those in the best model 
A: Model G with ?7l smaller by 10 km s~^ than that of the 
best model, model H with Ul larger by 10 km s^^ than that 
of the best model, model 1 with C/l larger by 2 km s~^ than 
that of the best model, and model J with Us larger by 10 
km s~^ than that of the best model. To investigate the dy- 



namical effects of the SMC on the orbital evolution of the 
LMC, wc investigate the model K without the SMC. 

Figure 2 describes the past ~ 9 Gyr orbital evolution of 
the Clouds in the best orbital model A. Here negative values 
of the time, T, represent the past, with T = corresponding 
to the present epoch. As shown in this figure, the present 
orbital period of the Clouds about the Galaxy is ~ 1.5 Gyr 
for the adopted gravitational potential and the masses of 
the Cloud and the Galaxy in the best model. Although the 
LMC-SMC distance remains very small (< 40kpc) over the 
last 4 Gyr (T > —4 Gyr), it cannot keep its binary status 
more than ~ 5 Gyr: Disintegration of the present-day binary 
orbit is inevitable in this model. 

Figure 3 describes the orbital evolution of the Clouds 
for the models C and D and thereby shows how the dy- 
namical friction between the Clouds due to the presence of 
the dark halo of the LMC influences the orbital evolution of 
the Clouds. The model D with Mlmc = 2.0 x 10^" Mq and 
without dynamical friction between the Clouds is exactly the 
same as the best model in GN. As shown in this figure, the 
Clouds can keep its binary status more than ~ 9 Gyr in the 
model D without dynamical friction of the Clouds whereas 
they cannot in the model C with dynamical friction. This 
suggests that the models with dynamical friction between 
the Clouds have difficulties in keeping the binary status of 
the Clouds for more than several Gyrs. In the Appendix 
A, we discuss the duration of the LMC/SMC binary status 
more extensively based on large number of orbital models. 
Figure 4 summarizes the orbital evolution for the represen- 
tative models, B, I, J, and F, for which we investigate the 
LMC evolution in the N-body simulations. 

2.2 N-body models 

2.2.1 The aelf-gravitating LMC 

The LMC is modeled as a fully self-gravitating system and 
composed of a live dark halo and a thin exponential disk 
with no bulge. The total mass of the dark halo, that of 
the disk, and the size of the disk are Mdm, and J?d, 
respectively. The mass ratio of the dark halo to the total 
mass (i.e., Mdm/MLMc) is fixed at 0.7 throughout the pa- 
per, which is consistent with the observation by vdMAHS. 
We adopt the Plummer potential in calculating the orbital 
evolution of the LMC (§2.1.1) and accordingly we need to 
adopt the corresponding density profile for the dark mat- 
ter halo of the LMC in this self-gravitating N-body models 
for self-consistency. The density profile of the dark halo is 
described as: 



(6) 



where tl is the distance from the center of the LMC and 
the value of ol is identical with that adopted in the equa- 
tion (2). Recent cosmological simulations within the frame- 
work of the CDM model (Navarro, Frenk & White 1996) 
have demonstrated the "universal" density distribution (the 
NFW profile); 

Po 



p(r) : 



(7) 



(r/rNFw)(l ± r/rNFw)^ ' 

where r, po, and rNPw are the spherical radius, the cen- 
tral density of a dark halo, and the scale length of the halo 
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Table 2. Model parameters for N-body simulations 
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respectively. For comparison, we also investigate the mod- 
els with this NFW profile for the LMC's dark matter halo, 
though it is not totally self- consistent with our orbital mod- 
els with the Plummer potential. 

The radial (7?) and vertical (Z) density profile of the 
initially thin disk of the LMC are assumed to be propor- 
tional to exp{—R/Ro) with scale length Ro = 2.6 kpc and 
to sech^(Z/Zo) with scale length Zo = 0.2J?o, respectively. 
We adopt the value of 2.6 kpc from Kinman et al. (1991) 
who investigated spatial distribution of RR Lyrae filed stars 
that are believed to be old, because our initial stellar disk 
is assumed to be older than ~ 9 Gyr. The circular velocity 
of the disk becomes a maximum value of Vm at = 5 kpc 
from the center of the disk for the adopted dark matter mass 
profile and Vm is 71 km s'^ . In addition to the rotational 
velocity made by the gravitational field of disk and halo com- 
ponent, the initial radial and azimuthal velocity dispersion 
are given to the disk component according to the epicyclic 
theory with Toomrc's parameter (Binncy & Trcmainc 1987) 
Q = 1.5. The vertical velocity dispersion at a given radius is 
set to be 0.5 times as large as the radial velocity dispersion 
at that point, as is consistent with the observed trend of 
the Galaxy (e.g., Wielen 1977). In order to compare phys- 
ical properties of now GCs formed from gas with those of 
old GCs initially within the LMC disk, we assume that the 
LMC disk initially contains 100 old GCs and the GC sys- 
tem has a disky distribution and rotational kinematics. We 
adopt the GC number of 100 that is much larger than the 
observed one (~ 13; van den Bcrgh 2000a), because we need 
an order of 100 GCs to evaluate structural and kinematical 
properties of the GC system: Only ~ 10 GCs do not allow 
us to derive the density profile and the rotational properties 
of the GC system. The assumption on the GC kinematics 
and structure is consistent with observations of structural 
and kinematics of old GCs in the LMC (e.g., Freeman et al. 
1983; Schommer et al. 1992). 

The disk is composed both of gas and stars with the gas 
mass fraction (/g) being a free parameter and the gas disk 



is represented by a collection of discrete gas clouds (corre- 
sponding to giant molecular clouds; GMCs) that follow the 
observed mass-size relationship (Larson 1981). Every pair of 
two overlapping gas clouds is made to collide with the same 
restitution coefficient of /r (Hausman & Roberts 1984). We 
vary the values of /r from 0.0 (no dissipation) 1.0 (highly 
dissipative) and thereby investigate the parameter depen- 
dences of the results on /r. We mainly present the results 
of the models with /r =0.5 and show some parameter de- 
pendences of the present results on /r. Although adopted 
method of "sticky particles" has been proven to be capable 
of addressing successfully some aspects of the hydrodynami- 
cal interactions in interstellar medium (ISM) for disk galax- 
ies (e.g., Hausman & Roberts 1984; Combes & Gerin 1985), 
it has some disadvantages in dealing with the more realis- 
tic physical processes of the ISM, such as hydrodynamical 
interaction between the hot interstellar gas (with the tem- 
perature of 10® K) and GMCs (e.g., evaporation of GMCs 
by the hot gas). We discuss this point in §3 for the results 
that may depend sensitively on the way to treat with ISM. 

In order to construct as a realistic gas disk model as 
possible, we consider the radial dependence of the gas mass 
fraction Fg{rL,) in the initial LMC disk. We expect that the 
inner gas mass fraction in the LMC is smaller than the outer 
one owing to more rapid consumption of gas in the inner 
regions with higher gas density. We therefore adopt the fol- 
lowing rule; 

Fg(rL) <x WC'^l) « « Sg«(rL), (8) 

Sg(rL) 

where tl, tsf, Sg, Eg, and a are the distance from the cen- 
ter of the LMC disk, the gas consumption time scale, the 
initial gas density, the gas consumption rate, and the pa- 
rameter controlling the radial dependence. Since we adopt 
the Schmidt law with the exponent of 1.5 for star formation 
(described below), the reasonable value of a is —0.5. Accord- 
ing to the value of Fg(rL,) derived from the above equation 
(8) at each radius, we determine the reasonable number of 
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old stellar particles and gaseous ones at each radius and 
thereby allocate these particles to each radial bin. By as- 
suming that the disk is composed only of gas initially, we 
determine Fgir-i) through the equation (8) and derive a rea- 
sonable radial distribution of gas and stars. To obtain more 
realistic initial stellar and gaseous distributions (e.g., with 
gaseous and stellar spiral arms), a LMC disk with Fg(rL) is 
allowed to relax for 10 dynamical time (~ 0.7 Gyr). We then 
use this disk as an initial LMC disk model in the simulations. 

2.2.2 Star formation and chemical evolution 

The gas is converted into either field stars or globular clus- 
ters (GCs), so that we distinguish the formation process of 
field stars from that of GCs throughout this paper. Field 
star formation is modeled by converting the coUisional gas 
particles into collisionless new stellar particles according to 
the algorithm of star formation described below. We adopt 
the Schmidt law (Schmidt 1959) with exponent 7 = 1.5 (1.0 
< 7 < 2.0, Kennicutt 1998) as the controlling parameter of 
the rate of star formation. The amount of gas consumed by 
star formation for each gas particle in each time step is given 
as: 

Pg oc pg^, (9) 

where pg is the gas density around each gas particle. We 
convert a gas particle into a field star only if the local gas 
density p exceeds the observed threshold gas density of ~ 3 
Mq pc~^ for Magellanic dwarf irregular galaxies (Hunter et 
al. 1998). These stars formed from gas are called "new stars" 
(or "young stars") whereas stars initially within a disk are 
called "old stars" throughout this paper. 

We use the cluster formation criteria derived by pre- 
vious analytical works (e.g., Kumai et al. 1993) and hydro- 
dynamical simulations with variously different parameters of 
cloud-cloud collisions on a 1-lOOpc scale (Bekki et al. 2004a) 
in order to model GC formation. A gas particle is converted 
into a cluster if it collides with other high velocity gas (with 
the relative velocities ranging from 30kms~^ to lOOkms"^) 
and having an impact parameter (normalized to the cloud 
radius) less than 0.25. Although both binary cluster and 
single one are formed during high-velocity cloud-cloud colli- 
sions (Bekki et al. 2004a), we assume that only one cluster is 
formed from one event of cloud-cloud collision for simplicity. 
These GCs formed from gas are called "new GCs" whereas 
GCs initially within a disk are called "old GCs" throughout 
this paper. 

This model is strongly supported by recent observations 

(e.g., Zhang et al. 2001) that have revealed that there is a 
tendency for young clusters to be found in gaseous regions 
with higher velocity dispersion, where cloud-cloud collisions 
are highly likely. In the present model, the clusters are not 
formed in the isolated model (described later) at all. This 
enables us to investigate how and whether the tidal inter- 
action between the Clouds and the Galaxy triggers the for- 
mation of GCs. About an order of 10^ clusters are formed 
in the models where the LMC interacts with the SMC and 
the Galaxy. The preset simulations only investigate dynam- 
ics with the scale down to ~ 100 pc so that it does not allow 
us to investigate which new clusters with the size of ~ 10 pc 
survive from tidal force of the Galaxy and the Clouds to be 
observed as GCs at the present time. Therefore, we assume 



that all clusters formed in the simulations become GCs, and 
thereby analyze the physical properties of new GCs. The to- 
tal number of GCs in the simulations may be overestimated 
in the present model because of this assumption. 

Chemical enrichment through star formation and su- 
pernovae feedback during the LMC evolution is assumed 
to proceed both locally and instantaneously in the present 
study. We assign the metallicity of original gas particle to 
the new stellar particle and increase the metals of the each 
neighbor gas particle with the total number of neighbor gas 
particles equal to A'gas, according to the following equation 
about the chemical enrichment: 

AMz = {ZiR^^ms+{l.Q-R^^t){l.Q-Zi)m^y^^t}/N^^{lQ) 

where the AMz represents the increase of metal for each gas 
particle. Zi, i?met, rua, and j/met in the above equation rep- 
resent the metallicity of the new stellar particle (or that of 
original gas particle), the fraction of gas returned to inter- 
stellar medium, the mass of the new star, and the chemical 
yield, respectively. The values of -Rmct , J/mct , and the initial 
metallicity are set to be 0.3 and 0.005, and 0.002, respec- 
tively. For these values, the final mean metallicity of the 
inner region of a LMC disk is consistent with the observed 
one with [Fe/H] = -0.3 ± 0.04 (Luck et al. 1998). 

In order to discuss the metallicity distribution of old 
stars stripped from the LMC disk, we need to assume that 
the old stellar disk has a metallicity gradient consistent with 
observations. Friel (1995) has derived the metallicity gradi- 
ent of the Galactic stellar disk based on the ages and metal- 
licities that are estimated for the Galactic open clusters. 
Since we do not have any available data on metallicity gra- 
dient of open clusters in the LMC, we compromise to use the 
observed slope of the Galactic metallicity gradient by Friel 
(1995) for the LMC old stellar disk. We therefore allocate 
metallicity to each disc star according to its initial position: 
at tl = Rl, where tl (Rl) is the projected distance in units 
of kpc from the center of the LMC disk, the metallicity of 
the star is given as: 

[m/H]^^^j,^ = [m/H]^_^^^o + «d x Rl. (11) 

If we adopt a plausible values of —0.091 for the slope aa 
(Friel 1995) and the central value of —0.73 for [m/H]^ rL=0' 
the mean metallicity of the LMC old disk is —1.0 in [Fe/H]. 

2.3 N-body evolution on the pre-determined 
orbits 

We numerically investigate the evolution of the LMC disk 

under the gravitational influence by the Galaxy and the 
SMC for the last ~ 9 Gyr (128 dynamical time scales of the 
LMC) by using the above N-body models of the LMC. We 
adopt an orbital model (e.g.. A) and consider that the LMC 
(the SMC) in a simulation is always on the predetermined 
orbit. Based on the orbits of the Clouds derived in the above 
orbital calculations (i.e., models A — K), we create a look- 
up table of positions and velocities of the Clouds at each 
time step for —8.8 < T < Gyr with the time step width 
of 1.4 X 10® yr in each orbital model. The center of mass 
in the LMC (the SMC) at each time step in a simulation 
with an orbital model is set to be the same as the location 
of the LMC (the SMC) at the time step in the look-up table 
of the orbital model. A fully self-gravitating LMC model is 
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Figure 5. Morphological evolution of old stars and gas seen from the face-on view (referred to as the x-y plane) with respect to the 
LMC disk for the fiducial model. 



influenced by the same fixed potential of the Galaxy and the 
SMC used in the equation (1) and (3), respectively, in a sim- 
ulation. Since we do not intend to constitute the SMC as a 
sclf-gravitating particle system in the present study, possible 
important physical processes such as direct hydrodynamical 
interaction between gaseous components of the Clouds and 
mass-transfer from the SMC to the LMC are not included. 

The initial spin of a LMC disk in a model is specified by 
two angles, 9 and <j}, where 9 is the angle between the Z-axis 
and the vector of the angular momentum of a disk and cj) is 
the azimuthal angle measured from X-axis to the projection 
of the angular momentum vector of a disk onto the X — Y 



plane. We mainly describe the results of the models with 9 
= 99° and cf) = 257°, because the final structural proper- 
ties of the LMC disk at T = Gyr in this model is broadly 
consistent with observations in that (1) the simulated LMC 
disk has a stellar bar with the size of the bar similar to the 
observed one, (2) the position angle of the simulated bar 
projected onto the sky is not largely different from the ob- 
served one (yet not exactly the same), and (3) the northeast 
side of the simulated disk is the near side with respect to the 
Galaxy, which is consistent with the observations by Cald- 
well & Coulson (1986). The LMC disk precesses and nutates 
under the tidal torque from the Galaxy during the dynam- 
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Figure 6. The same as Figure 5 but for new stars formed from gas. 



ical evolution of the LMC (Weinberg 2000). Therefore it is 
very hard to choose the initial values of 6 and 4> for which the 
final internal spin axis of the LMC at T = in our simula- 
tions is nearly the same as those inferred from observations. 
Thus, we show the results of the models for which the fi- 
nal structural properties of the LMC are broadly consistent 
with observations. 

We mainly describe the "fiducial model" with the best 
orbital model/typo k, 6 = 99°, 6 = 257°, /g = 0.5, and 
= 0.5. This is firstly because the parameters of this model 
are the most reasonably consistent with observations, and 
secondly because this model shows typical behaviors in dy- 



namical and chemical evolution of the LMC influenced by 
the Galaxy and the SMC. We also show the results of the 
representative 21 models thereby discuss the dependences of 
the results on model parameters such as the orbital types, 
9, 4>, /g, and /r. The values of these parameters are summa- 
rized for each model in the Table 2: Model number (column 
1), orbital type (2), in units of degrees (3), <f) in units of 
degrees (4), and the comments on the models (5). 

The initial particle number used in a self-gravitating 
LMC model is 50000 for the dark matter, 25000 for the old 
stars, and 25000 for the gas. The total particle number is 
increased up to ~ 160000 owing to the formation of new 
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field stars and globular clusters. All the simulations have 
been carried out on GRAPE board (Sugirnoto ot al. 1990) 
at the Astrononiical Data Analysis Center (ADAC) of the 
National Astronomical Observatory of Japan. The parame- 
ter of gravitational softening is sot to be fixed at 0.15 kpc. In 
the following, in order to show more clearly the morphologi- 
cal and kinematical properties of the simulated LMC, we set 
the face-one view (edge-on view) of the LMC disk to be al- 
ways the x-y (x-z) plane by rotating the LMC disk by some 
degrees at a given time T. Thus it should be noted in the 
following that the x-y {x-z and y-z) plane is not identical 
with the X-Y {X-Z and Y-Z, respectively) in the Figure 1. 



3 RESULTS 

3.1 The fiducial model 

3.1.1 Formation of bar and thick disk 

Figures 5, 6, and 7 describe the morphological evolution of 
the LMC for the last ^ 9 Gyr in the fiducial model. Owing to 
the smaller mass of the LMC (Mlmc = 10^° Mq), the LMC 
disk is influenced by the strong tidal field of the Galaxy from 
the early dynamical evolution of the LMC. During the first 
pericenter passage of the LMC with respect to the Galaxy (T 
~ —6.8 Gyr), two spiral arms composed of old stars and gas 
are formed within the disk owing to the tidal perturbation 
from the Galaxy. The distance between the Clouds is larger 
than ~ 100 kpc for T < -6 Gyr so that the SMC does 
not dynamically influence the LMC in the early dynamical 
evolution of the LMC. After the first pericenter passage, the 
bar-like structure composed mostly of old stars is formed in 
the central region of the disk (T = —5.5 Gyr). 

Formation of field stars starts gradually from the cen- 
tral region of the disk, in particular, from the high density 
regions of gas along the inner spiral arms of the disk (T 
= —7.1 Gyr). As the LMC first passes by the pericenter 
of its orbit with respect to the Galaxy at T = —6.8 Gyr, 
the gaseous regions where field stars are actively forming 
are shifted from the center to the two remarkable spirals 
formed by the Galactic tidal perturbation. After the forma- 
tion of the central bar- like structure {T = —5.5 Gyr), the gas 
density within the bar-like structure becomes higher owing 
to the enhanced rates of cloud-cloud collisions along/within 
the bar-like structure. As a result of this, the field star for- 
mation becomes more efficient within the central bar-like 
structure. The LMC disk is thus morphologically classified 
as a barred spiral in the early phase of its evolution (T < 
—5.5 Gyr), where the SMC does not dynamically influence 
the LMC. 

As the strong tidal interaction between the Clouds 
starts (T = —3.8 Gyr), the stellar bar grows and thus be- 
comes more remarkable compared with the outer stellar disk 
(T = —3.3 Gyr). The formation of the stronger bar is asso- 
ciated closely with (1) the tidal perturbation from the SMC 
and (2) formation of new stars along/within the bar. A sig- 
niflcant fraction of old stars (~ 17 %) initially within the 
outer part of the disk are tidally stripped and field star for- 
mation is still ongoing mostly within the bar after the first 
LMC-SMC encounter. Consequently, most of old and new 
stars are located in the central bar of the disk at T = —2.2 
Gyr. Tidal interaction between the Clouds and the Galaxy 
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Figure 7. The same as Figure 5 but seen from the edge-on (the 
X-Z plane) for four different epochs. For comparison, the final 
mass distributions of old stars (left) and new stars (right) seen 
from the cdgc-on view in the isolated disk model at T = Gyr 
are shown in the bottom two panels. Note that owing to the tidal 
stripping of outer halo stars by the Galactic strong tidal field, the 
LMC slowly develops its stellar halo with an inhomogeneous den- 
sity distribution. Note also that after the stellar bar formation (T 
= —5.5 Gyr), a very thick disk is finally formed. The thick disk 
is significanly thicker than the old stellar disk and the new stellar 
one in the isolated model, which implies that tidal intera<;tion be- 
tween the Clouds and the Galaxy is a main cause of the thick disk 
formation (Compare the fiducial model with the isolated one). 



finally results in the formation of an elliptic disk which sur- 
rounds the central bar, consists mainly of old stars, and has 
a major axis misaligned significantly with that of the central 
bar (T = Gyr). 

As is shown in Figure 7, the initially thin stellar disk 
is finally transformed into a thick one with the size signif- 
icantly smaller than its original one owing to the tidal in- 
teraction between the Clouds and the Galaxy for the last ~ 
9 Gyr. The strong Galactic tidal field dynamically heats up 
the stellar disk every time the LMC passes by the pericenter 
of its orbit with respect to the Galaxy (e.g., T ~ —6.8 Gyr). 
The vertical heating by the Galaxy and the resultant thick- 
ening of the LMC disk starts from the outer, more fragile 
part of the disk. After the stellar bar formation (T = —5.5 
Gyr) , the disk becomes severely warped and shows a sign of 
a "banana" shape at T = —3.8 Gyr. The bending instability 
following bar formation in tidal galaxy interaction also con- 
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Figure 8. The time evolution of the lialf-mass radius Rh and the 
vertical velocity dispersion (ctz) for the fiducial model (solid) and 
the isolated one (dotted). The value of R^^ (ctz) at each time step, 
denoted as Rh{T) {az{T)), is normalized to the initial value of 
Rh{0) (o-z(O)) at T = -8.8 Gyr. Note that Rh in the fiducial model 
dramatically decreases owing to the central new stars efiiciently 
formed during the tidal interaction between the Clouds and the 
Galaxy. Note also that (1) CTz also significantly increases for —6 
< T < —2 Gyr in the fiducial model and (2) the degree of the 
increase is much more significant in the fiducial model than in 
the isolated one. These two results suggest that the increase of 
(Tz is not due to the numerical heating caused by small number 
of stellar particles but due to the tidal heating from the Galaxy 
and the SMC. 
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tributes significantly to the thickening of a disk (e.g., Miwa 
& Noguchi 1999). The thickening of the LMC disk seen in 
Figure 7 for —3.8 < T < Gyr is due partly to the bending 
instability in the disk. 

As the disk thickening continues, old stars stripped from 
the outer part of the disk (i?L > 5 kpc) gradually develops 
the stellar haJo around the LMC disk. Although the distri- 
bution of stars in the halo is quite elongated and inhomo- 
geneous at T = —3.8 Gyr, it becomes more spherical and 
homogeneous at T = Gyr. About 17 % of the initial old 
stars are spatially redistributed to form the outer stellar halo 
and only 2 % of the halo is composed of new stars formed 
from gas and then tidally stripped during the interaction 
between the Clouds and the Galaxy. One of the important 
characteristics of the stellar halo is that the mass fraction 
of the halo to the disk finally becomes an order of ~ 10 %, 
which is significantly higher than that of the Galaxy (a few 
percent; e.g.. Freeman 1987). The physical properties of the 
stellar halo formed from the Galaxy-LMC-SMC interaction 
will be discussed later. 

Figure 8 summarizes the time evolution of global struc- 
tural and kinematical properties of the LMC stellar disk for 
the last ~ 9 Gyr. Firstly, the effective radius represented by 
Rh (the half-mass radius) of the disk becomes significantly 
small owing to the formation of new stars in the central re- 
gion of the disk. For example, Rh shows 65 % of the orig- 
inal value for T = —3.8 Gyr (before the commencement of 
the strong tidal interaction with the SMC) and finally shows 
^ 33 % of the original one at T = Gyr: The disk size be- 
comes smaller as the disk is dominated more significantly by 
young stellar populations. The isolated model shows only ~ 
22 % decrease of its disk size owing to the much less efficient 
star formation in its central region, which confirms that the 
dramatic change of the size of the LMC disk in the fiducial 
model is caused by central star formation triggered by the 
Galaxy-LMC-SMC interaction. 

As is shown in Figure 8, the vortical velocity dispersion 
represented by tJz dramatically (a factor of ^ 3) increases 
during the Galaxy-LMC-SMC interaction for the last ~ 9 
Gyr, whereas a^, of the disk in the isolated disk model only 
gradually increases by only a factor of ~ 1.4. The increase in 
CTz for the isolated model is caused by (1) the scattering of 
stars via gas clouds, (2) "self-heating" by the spiral arms in 
the disk (e.g., Sellwood & Carlberg 1984), and (3) numerical 
heating due to the adopted small particle number (an order 
of 10^) in the present simulation. Given the fact that the 
increase in ctz in the fiducial model is more than a factor of 
2 larger than that in the isolated one, the contribution of 
numerical heating to ctz in the fiducial model is quite minor. 
Therefore, it is reasonable to say that the increase in ctz in 
the fiducial model is caused mainly by the tidal heating from 
the Galaxy and the SMC. 



Time (Gyr) 

Figure 9. Time evolution of cloud-cloud collision number (up- 
per), field star formation rate (middle), and globular cluster for- 
mation rate (lower) in the fiducial model. 



3.1.2 Stellar populations 

Figure 9 shows the formation histories of field stars and 
GCs in the fiducial model. The formation rate of field stars 
gradually increases to reach ~ 0.1 Mq yr^^ in the first ~ 1 
Gyr dynamical evolution owing to the development of high 
density gas regions along the spiral arms in the disk. After 
the first pericenter passage of the LMC (T ~ —6.8 Gyr), 
the star formation rate is moderately enhanced and reaches 
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Figure 10. The radial gradients of age (lower) and metallicity 
(upper) in the LMC disk at T = Gyr in the fiducial model. 



T= Gyr 




^ 0.4 Mq yr~^ (averaged for 0.1 dynamical time of the 
LMC corresponding to ~ 13 Myr) owing to the formation 
of strong, double-armed gaseous arms, where gas density 
becomes significantly high. The cloud-cloud collision rate is 
dramatically (more than an order of magnitude) enhanced 
within ~ 0.5 Gyr after the perieenter passage, because the 
strong Galactic tidal effects increase the velocity dispersion 
of the gas clouds. However, the cloud-cloud collisions that 
leads to the formation of GCs in the disk do not occur in this 
first perieenter passage, just because the tidal perturbation 
is not strong enough to trigger cloud-cloud collisions with 
moderately high relative velocities (between 30 and 100 km 
s~^) and small impact parameters (< 0.25). 

Owing to the rapid gas consumption by formation of 
field stars during/after the first perieenter passage of the 
LMC, the formation rate of the field stars does not increase 
significantly until the strong Galaxy-LMC-SMC tidal inter- 
action begins at T ~ —3.8 Gyr. The tidal perturbation from 
the SMC and the Galaxy triggers the moderately enhanced 
star formation rate of ~ 0.1 Mq for —3.5 < T < —2 Gyr. 
During this period, the GC formation also becomes efficient, 
essentially because the combined tidal efi^ects of the Galaxy 
and the SMC arc strong enough to enhance cloud-cloud col- 
lisions required for GC formation in the LMC disk. The peak 
of the GC formation (~ T ~ —2.5 Gyr) is nearly coincident 
with that of the field stars for -3.5 < T < -2 Gyr. The GC 
formation also occurs at T ~ —0.2 Gyr, when the LMC-SMC 
distance becomes very small (less than 10 kpc; smaller than 
the original LMC disk size) so that the LMC collides with 
the SMC. This final collision between the Clouds around 
0.2 Gyr ago could have significant effects in the recent star 



Figure 11. Distribution of all particles seen from the face-on 
view of the LMC disk (upper) and the radial distribution of the 

mass fraction of new stars to all stars (lower) at T = Gyr in the 
fiducial model. The mass fraction of new stars is estimated along 
the major axis of the stellar bar (solid) and along the minor axis 
(dotted). 



formation histories of the LMC, as suggested by previous 
authors (e.g., GSF). 

About 47 % of the initial gas is converted into new field 
stars for the last ~ 9 Gyr evolution of the LMC and most of 
new stars are concentrated in the central bar. Only 0.5 % of 
the gas is converted into GCs, which reffects the fact that the 
cloud-cloud collisions with moderately high speed (between 
30 and 100 km s~^) and small impact parameter (< 0.25), 
required for cluster formation, do not occur until the LMC 
begins to interact violently with the SMC when the two are 
less than 10 kpc apart (T=— 3.6Gyr). Chemical enrichment 
resulting from the moderately enhanced star formation in- 
creases gradually the metallicity of new field stars and GCs. 
For the fiducial model with the initial gaseous metallicity of 
0.002 ([Fe/H] = -1.0), j/met (chemical yield) of 0.005, and 
Rmet of 0.3, the mean metallicity of new field stars within 
the disk and the halo finally becomes 0.007 (-0.46 in [Fe/H]) 
at T = Gyr (It should be emphasized here that the final 
metallicity depends strongly on the initial values of the ini- 
tial metallicity, j/met, and Rmet). 

Chemical enrichment proceeds more in the central re- 
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Figure 12. The age distribution of new stars inside the stellar 
bar (solid) and outside the bar (dotted) at T = Gyr in the 
fiducial model. A new star is regarded as "inside the bar" , if the 
distance of the star along the major axis of the bar is less than 
5 kpc and if that along the minor axis are is less than 1.5 kpc. 
For convenience, the normalized number of new stars is shown. 
The epochs of the pcriccntcr passage of the SMC with respect to 
the LMC are shown by thick arrows for comparison. Note that 
the age distribution of new stars inside the bar shows the peak 
around ^ 2 Gyr whereas that outside the bar shows the peak 
around 6 Gyr. 



gion where gas consumption by star formation is more rapid 
and efficient owing to the gas inflow triggered by the Galaxy- 
LMC-SMC interaction. The disk consequently shows a neg- 
ative motaUicity gradient for new field stars in the sense that 
the inner regions show higher mctallicity than the outer ones 
at T = Gyr (Figure 10). The disk also shows a positive 
age gradient with the new field stars within the central 1 
kpc of the disk being a few Gyr younger than those outside 
the central 1 kpc. Both age and mctallicity gradients are 
rather flat for 3 < -Rl < 7 kpc in the disk (where f?L is the 
projected distance from the center of the LMC disk), pos- 
sibly because the stellar bar dynamically mixes the young 
stellar population with different ages and metallicities owing 
to the stream motion of stars and gas along the bar. Both 
the age and metallicity dispersions are larger in the cen- 
tral regions, which reflect the fact that gaseous components 
with different metallicities are transferred to the center and 
converted to new stars there in different epochs during the 
Galaxy-LMC-SMC interaction. It should be stressed here 
that recent observations have found that intermediate-age 
stars are more centrally concentrated than the older stars in 
the LMC disk (e.g., Cole et al. 2000). 

Figure 11 clearly shows that the mass fraction of new 
field stars (referred to as fns hereafter) along the major axis 
of the bar is significantly larger within the stellar bar than 
outside it. About 50 % of the stars within the central ~ 500 
pc of the disk along the major axis of the bar are new field 
stars formed from gas during the Galaxy-LMC-SMC interac- 
tion. The radial profile of fns shows a relatively flat density 
profile ("shoulder") with /ns ~ 0.3 (0.2) along the major 
(minor) axis of the bar and has abrupt steepening of the 
gradients around |_Rl| ~ 1 kpc and 4 kpc. The radial gradi- 
ent for |i?L| < 7 kpc appears to be slightly steeper along the 
minor axis of the bar than along the major one. /ns becomes 
smaller than 0.1 for 5 kpc < |ilL|, which indicates that the 




-5 
X-axis (kpc) 

Figure 13. Radial dependences of line-of-sight- velocity (Vy) and 
velocity dispersion (cry ) along with the x-axis for old stars (upper) 
and for new stars (lower) at T = Gyr in the fiducial model. 
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Figure 14. Radial dependences of vertical velocity dispersion 
((Tz) for old stars (thick solid) and new stars (thin solid) with the 
vertical distance \z\ from the disk plane of the LMC less than 1 
kpc (upper) and those with \z\ less than 7.5 kpc (lower) at T = 
Gyr in the fiducial model. For comparison, the results of the 
isolated model are also shown by dotted lines. 



outer disk is dominated by old stellar populations. These 
radial dependences of fns are due essentially to the centrally 
concentrated young stellar populations formed during the 
Galaxy-LMC-SMC interaction. 

Figure 12 shows that the stellar age distribution is sig- 
niflcantly different between new stars inside the bar and 
those outside the bar. The stellar age distribution inside the 
bar has a peak around ~ 2 Gyr, which means that the stel- 
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Figure 15. The distribution of new stars on the age- 14 plane, 
where 14 is the line-of-sight velocity parallel to the z-axis (upper) 
and the mean value of the vertical velocity dispersion (Tz as a func- 
tion of age for the new stars (lower) at T = Gyr in the fiducial 
model. The lower figure therefore describes the age-dispersion re- 
lation of new stars. For comparison, the age-dispersion relation 
for the isolated model is shown by a dotted line in the lower panel. 

lar population inside the bar is dominated by new stars that 
are formed after the Galaxy-LMC-SMC interaction becomes 
stronger owing to the dynamical coupling of the Clouds. The 
bump around the ages of 6 — 8 Gyr in the age distribution in- 
side the bar suggests that some fraction of new staxs formed 
in the early Galaxy-LMC interaction. The age distribution 
outside the bar, on the other hand, has a peak around 6 Gyr, 
which means that the stellar population outside the bar is 
dominated by new stars formed in the early evolution of the 
LMC, in particular, those formed during the Galaxy-LMC 
interaction 6 — 7 Gyr ago. These results in Figures 11 and 
12 suggest that (1) the dominant stellar population is signif- 
icantly different between different regions of the LMC disk 
and (2) the difference is due essentially to the Galaxy-LMC- 
SMC interaction which forms a bar and thus drives efficient 
inner transfer of interstellar gas. 

3.1.3 Kinematics 

The long-term tidal perturbation from the Galaxy and the 
SMC to the LMC causes dramatic changes in kinematical 
properties of the disk that is initially "dynamically cold". 
Figure 13 shows that the rotational velocity of old stars 
(hereafter referred to as Vrot, and represented as Vy in Fig- 
ure 13 for convenience) has a peak value of ~ 50 km s~^ 
around the central 3 — 4 kpc and then decreases sharply to- 
ward outward from there. The ratio of the maximum value 
of Kot to Vm (the maximum value of the circular velocity Vc 
in the initial disk) for old stars is ~ 0.7, which implies that 
azimuthal/radial velocity dispersion of old stars is signifi- 
cantly increased by the Galaxy-LMC-SMC interaction. Vrot 
decreases outside the central 3 — 4 kpc much more sharply 
than expected from the initial mass profile of the LMC. The 



initial disk shows a decrease by a factor of ~ 17 % in 14: 
(and thus Vrot) for 5 < Rh < 10 kpc whereas Viot in the fi- 
nal LMC disk shows a decrease by a factor of ~ 60 % in Kot 
for the corresponding region. This clearly indicates that the 
long-term dynamical heating of the LMC disk by the Galaxy 
and the SMC changes the shape of the rotation curve of the 
LMC. 

The initial radial profile of azimuthal velocity dispersion 
of old stars is monotonously decreasing toward the outer 
region of the LMC disk. As shown in Figure 13, the velocity 
dispersion ay (a measure of azimuthal velocity dispersion) 
shows a large value of 30 — 40 km s~^, which is even larger 
than the central one of ~ 20 km s~^ , around 7 — 8 kpc from 
the center of the disk. This result confirms that the long- 
term dynamical heating by the Galaxy and the SMC drives 
the kinematical change of the LMC's outer disk that is more 
susceptible to external tidal perturbation. The results of the 
new stars are essentially the same as those of old ones, except 
that (1) the velocity dispersion ay is on average smaller in 
new stars than in old ones and (2) the maximum Viot is 
only slightly larger in new stars than in old stars probably 
because of gaseous dissipation. The asymmetric profiles of 
ay and Kot seen both in old stars and in new ones are due 
partly to the collision between the Clouds at T ^ —0.2 Gyr. 

Figure 14 shows that the radial gradient of the vertical 
velocity dispersion cTz of the disk for stars with the vertical 
distance {\z\) from the disk plane equal to or less than 1.5 
kpc (i.e., those stars within the thick disk). It is clear from 
this Figure 14 that a^ of old stars shows a monotonous de- 
crease for J?L < 4 kpc and begins to increase gradually and 
slightly toward outward from there, though the profile for 
-Rl > 4 kpc is somehow irregular. The outwardly increasing 
(Tz is more remarkable for old stars with \z\ < 7.5 kpc, be- 
cause these stars include some fraction of old halo stars that 
form a kinematically hot stellar system around the LMC. 
The larger a^ in the outer LMC disk means that the scale 
height also increases with radius for _Rl > 4 kpc. These re- 
sults arc broadly consistent with recent observational results 
by Alves & Nelson (2000) and vdMAHS. 

fTz is increased by a factor of 2.4 in the center of the disk 
and a factor of 5.5 at _Rl = 8 kpc for the old stars with \z\ 
< 1.5 kpc owing to the long-term tidal interaction between 
the Galaxy and the Clouds. The radial profiles of ctz for new 
stars also show such an outwardly increasing for Rl, > 4 
kpc, though the increase is less remarkable compared with 
that seen in old stars. CTz of new stars is on average lower than 
that of old stars for nearly every radii, essentially because 
new stars are formed from gas that dissipates away random 
kinematical energy yielded by the Galaxy-LMC-SMC inter- 
action. The central "dip" in the radial profile of a^ for new 
stars is also caused by efficient gaseous dissipation there. 
Thus Figures 13 and 14 suggest that the outer kinematics of 
the LMC disk have valuable information on the past inter- 
action history of the LMC with the Galaxy and the SMC. 

Figure 15 describes a relation between velocity disper- 
sion of CTz and ages (i.e., so-called "age-dispersion relation") 
for new stars of the LMC disk. The isolated model shows the 
trend of decreasing a^ with decreasing ages, which refiects 
the fact that a new star formed more recently originates 
from gas that has experienced a larger amount of gaseous 
dissipation. This tendency is seen in the disk of the fiducial 
model for the new stars with the ages older than ~ 4Gyr, 
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Figure 16. The distribution of halo old stars (upper left) and 
halo new ones (lower left) seen from the edge-on of the LMC disk 
and the radial dependences of the linc-of-sight velocity parallel to 
the y-ax\s (Vy for halo old stars (upper right) and for halo new 

stars (lower right) at T = Gyr in the fiducial model. Here a star 
can be classified as a "halo" star either if the projected distance 
from the center of the LMC disk is less than 7.5 kpc or if the 
vertical distance from the LMC disk plane is less than 1.5 kpc. 
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Figure 17. The projected density distribution of halo stars at 
T = Gyr (solid) and at T = -5 Gyr (dotted) in the fiducial 
model. For comparison, the radial distribution of stars in the ini- 
tial exponential disk is shown by a dashed line for comparison. 
The definition of "halo stars" is the same as that described in 
Figure 16. 



though the mean Uz is more than a factor of 2 higher in the 
fiducial model than in the isolated model owing to the dy- 
namical heating by the Galaxy. However, the new stars with 
the ages younger than ~ 4 Gyr show an interesting tendency 
of (Tz slightly increasing with decreasing ages. This is prob- 
ably because a new star formed more recently originates 
from gas that is more strongly randomized by the combined 
tidal effects of the Galaxy and the SMC: The net effect of 
gaseous dissipation is weaker than that of randomization of 
gaseous motion by the stronger tidal perturbation from the 
Galaxy and the SMC for the last ~ 4 Gyr. These results 
imply that the age-dispersion relation in stellar populations 
of the LMC disk is significantly different from that observed 
in the Galaxy. 
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Figure 18. Final radii of halo old stars as a function of the initial 
radii (upper) and the number distribution of the initial radii of 
the stars that are regarded as "halo stars" at T = Gyr (lower) 
in the fiducial model. The definition of "halo stars" here is the 
same as described in Figure 16. The long-dashed and short-dashed 
lines represent the initial tidal radius rt (= 13 kpc) and ~ 3Ro 
(7.5 kpc), where Rq is the scale radius of the initial LMC disk, 
respectively. Stars above the dotted line are those for which the 
final radii are larger than the initial ones (i.e., they are transfered 
outward due to the tidal effects of the Galaxy and the SMC). 
For convenience, the normalized number (i.e. imiiiber fraction) 
is shown in the lower panel. Note that most of the halo stars 
originate from the outer part {R > 4 kpc) of the initial LMC 
disk. 
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3.1.4 Stellar halo properties 

One of the important outcome of the long-term tidal inter- 
action between the Clouds and the Galaxy is the formation 
of a stellar haJo through the redistribution of old stars ini- 
tially within the LMC disk. It should also be stressed here 
that our disk models initially do not include very old (~ 
10 — 13 Gyr), metal-poor ([Fe/H] < —1.6) stellar halo com- 
ponents that the Galaxy is observed to have. Therefore, the 
halo properties in the present study are dependent totally 
on those of old stars and new stars stripped from the original 
disk. Accordingly, direct comparison of the present numeri- 
cal results with observations is not simply possible, because 
the LMC stellar halo may contain the old stellar halo that 
was formed prior to its disk formation. Wc thus suggest that 
the present results may well be able to be compared with 
observational results for halo stars with ages less than ~ 9 
Gyr. 

As have been shown in Figure 7, some fraction of old 
stars stripped from the outer part of the disk due to the 
Galaxy-LMC-SMC interaction are redistributed in the outer 
halo region of the disk. Figure 16 shows structural and kine- 
matical properties of the stellar halo composed of stars ei- 
ther outside the disk (i.e., i?L > 7.5 kpc) or well above/below 
the thick disk (i.e., \z\ > 1.5 kpc). Although these old halo 
stars have a relatively homogeneous distribution within the 
central ~ 6 kpc, they have an elongated distribution for a 
wider field (12 kpc x 12 kpc) of view owing to the old stars 
being now dynamically influenced by the tidal fields of the 
Galaxy and the SMC. The new stars in the halo also show 
such an elongated distribution, though it is less remarkable 
compared with old halo stars. The new halo stars appear 
to be more flattened compared with old halo ones for the 
central ~ 6 kpc because of the smaller number of the stars 
with \z\ > 4.0 kpc. Since the mass fraction of the new stars 
in the halo region is only 2 %, the entire distribution of 
the halo is determined by old halo stars: The halo has an in- 
ner homogeneous distribution and an outer elongated (thus 
inhomogeneous) one. 

As shown in Figure 16, both old and new halo stars 
show a sign of rotation, in particular, for those within the 
central 2 kpc, if they are seen from the edge-on view. Al- 
though the inner stellar components with a certain amount 
of rotation should be regarded as parts of the thick disk 
rather than the halo, the moderately rotating inner halo is 
regarded as an important characteristic that the stellar halo 
has if it originates from the outer part of the initially thin 
stellar disk through the redistribution of the stars during the 
Galaxy-LMC-SMC interaction. Figure 17 shows that (1) the 
projected density of the stellar halo at T = Gyr is roughly 
approximated as an exponential proflle with the slope shal- 
lower than the original exponential disk for Rj, > 2 kpc, (2) 
the projected density of the stellar halo is on average more 
than an order of magnitude lower than the original stellar 
disk and the difference in the density depends strongly on 
the radius, and (3) the projected density profile changes sig- 
nificantly with time at a given radius and it becomes steeper 
as the time passes by. 

Figure 18 demonstrates that about 75 % of old stars 
in the halo (at T = Gyr) originate from the disk regions 
with the distances larger than 5 kpc. This is because the 
tidal stripping of old disk stars and the subsequent redistri- 
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Figure 19. The metallicity distribution function (MDF) of halo 
old stars is shown by a solid line for T = Gyr in the fiducial 
model. For comparison, the initial MDF of the LMC disk is shown 
by a dotted line. 

bution of the stripped stars more efficiently occur during the 
Galaxy-LMC-SMC interaction. As a natural result of this, 

the final metallicity distribution function (MDF) in the old 
halo stars is expected to be significantly different from that 
of the initial disk. Figure 19 shows a possible MDF of the 
old halo stars for the fiducial model with an initial nega- 
tive metallicity gradient (described in the equation (11) in 
§2) and the mean metallicity of [Fe/H] = -1. The MDF in 
this Figure 19 clearly shows a peak around [Fe/H] = —1.2 
~ —1.3, which is by ~ —0.3 (in dex) lower than that of the 
initial disk and is by ~ —0.4 (in dex) higher than that of 
the observed Galactic stellar halo (e.g.. Freeman 1987). The 
derived lower metallicity of the old halo stars is due to the 
stars originating preferentially from the outer disk regions 
where the stellar metallicity is lower owing to the negative 
metallicity gradient. 

The derived difference in the MDF between the sim- 
ulated stellar halo and the observed Galactic stellar halo 
probably reflects the fact that the formation history of the 
Galactic stellar halo is significantly different from what is 
described in this paper (e.g., Bekki & Chiba 2000, 2001). 
One of the significant differences between the simulated stel- 
lar halo and the observed Galactic one is that the simulated 
halo includes a fraction of relatively metal- rich (— 1 < [Fe/H] 
< —0.3) and moderately young stars, though the mass frac- 
tion of these populations among the entire halo population is 
very small (~ 2 %). These metal-rich halo components result 
from the tidal stripping of new stars formed the LMC disk in 
the later phase of the Galaxy-LMC-SMC interaction. There- 
fore, the detection of such metal-rich, young stars within the 
LMC halo region in future observations could be an evidence 
that the Galaxy-LMC-SMC interaction is partly responsible 
for, at least, some part of the LMC stellar halo. 

3.1.5 Properties of GCs 

Figure 20 summarizes structural and kinematical properties 
of old and new GCs in the fiducial model at T = Gyr. 
As shown in this Figure 20, the final spatial distribution is 
remarkably different between the old GCs and the new ones. 
The new GCs has a more compact distribution than the old 
GCs and most of them are concentrated within the central 
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Figure 20. The distribution of old GCs projected onto the x-y 
plane (top left) and onto the x-z one (top right) and tliat of new 
GCs projected onto the x-y plane (middle left) and onto the x-z 
one (middle right) at 7' = Gyr in the fiducial model. The bottom 
panel shows the radial profiles of vertical velocity dispersion (uz) 
for old GCs (thin solid) and for new GCs (thick solid). 
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Figure 21. Age distributions of field stars (dotted) and clusters 
(solid) at T = Gyr in the fiducial model. For convenience, the 
normalized fraction of stars in each age bin is shown. 



^ 3 kpc. This is essentially because the new GCs form only 
in the high-density gaseous regions where the cloud-cloud 
collisions with moderately high speed (between 30 and 100 
km s~^) and small impact parameter (< 0.25), required for 
cluster formation, occur. The face-on distribution of the now 
GCs is fairly elongated with the direction of the elongation 
nearly parallel to the major axis of the stellar bar: Most of 
the new GCs axe within the stellar bar in the disk. Nearly all 
of the new GCs is confined within 1 kpc from the LMC disk 
plane, which reflects the fact that the new GCs originate 
from the high-density regions in the thin gaseous disk. 

The old GCs, on the other hand, show a more wide 
spread distribution and arc not necessarily confined in the 
central bar region of the disk. The edge-on distribution of 
the old GCs shows a thick-disk appearance with the two GCs 
having \z\ of > 3 kpc and being able to be regarded as halo 
GCs. These two halo GCs are initially in the outer stellar 
disk and thus tidally stripped and subsequently distributed 
in the halo region. Figure 20 also shows a few differences 
in the final radial profile of vertical velocity dispersion (ctz) 
between the two different GC populations, though the final 
profiles of (Tz are severely influenced by the tidal perturba- 
tion from the Galaxy and the SMC, in particular, from the 
coUision with the SMC around T ~ -0.2 Gyr (Here the 
center of mass of jield stars is adopted for the estimation 
of kinematics). Uz is systematically higher in the old GCs 
than in the new GCs throughout the disk, because the old 
GCs experienced for a much longer time (~ 9 Gyr) the tidal 
heating from the Gala^xy (and the SMC) than the new GCs 
all of which are formed relatively recently (T < —3.3 Gyr). 

Figure 21 describes the difference in the age distribu- 
tions between the new field stars and the new GCs in the 
disk. It is clear from this figure that all clusters have ages 
younger than ~ 3.3 Gyr whereas the field star population 
show a wide distribution of ages. This result reflects the 
fact that the field star formation is sensitive to local gas 
density whereas the cluster formation occurs only when ran- 
dom motion of gas in the LMC becomes significantly large. 
The mean metallicity of the new GCs at T = Gyr is by 
—0.08 dex (in [Fe/H]) higher than that of the new field stars, 
because chemical enrichment associated with field star for- 
mation proceeds efficiently prior to the formation of new 
GCs (i.e., T < ^ —3 Gyr) and consequently the new GCs 
are thus formed from more metal-rich gas. The derived dif- 
ference suggests that if future observations on the detailed 
age distribution of field stars for the entire disk region of 
the LMC reveal the differences in age and metallicity dis- 
tributions between field stars and GCs, the difference can 
be understood in terms of the difference in the formation 
processes between field stars and GCs during the Galaxy- 
LMC-SMC interaction. 



3.2 Dependences on model parameters 

We have described the results of the fiducial model in which 
the model parameters are consistent with the latest observar- 
tions and thus the most realistic and reasonable. However, 
there are still some observational uncertainties in estimat- 
ing the total mass of the LMC (vdMAHS). Considering this 
uncertainty, we mainly show the results of the more massive 
model (model 9) with larger LMC mass (Mlmc = 2.0 x 
10^° Mq) in this subsection. We also briefly summarize the 



© 1994 RAS, MNRAS 000, 1-?? 



18 K. Bekki, M. Chiba 



10: T=,-B.SGyr 



■10 

10 ; 'S-S 



•10 -5 

X-3K^^ (kpc} 



10 -10 -5 S 10 



Figure 22. The same as Figure 7 but for the model 9. 




-6 -4 -2 
Time (Gyr) 



Figure 23. The same as Figure 8 but for the model 9. 



parameter dependences on less important parameters (e.g.. 
orbits of the Clouds and /r). 



3.2.1 LMC mass 

The adopted LMC mgiss Mlmc is consistent with the obser- 
vation by vdMAHS but smaller than that adopted in most 
of previous studies with Mlmc = 2.0 x 10^° Mq (GN). 
Accordingly, it is important to investigate how the LMC 
evolution depend on its mass if it heis a larger mass. The 
stellar disk of the more massive LMC in the model 9 (10) is 
less susceptible to tidal perturbation from the Galaxy and 
the SMC than that in the fiducial model, because it is more 
strongly self-gravitating. As a result of this, the SMC in this 
model does not so strongly influence the LMC disk dynam- 
ically even after it dynamically couples with the LMC (T ~ 
—6 Gyr) in comparison with the fiducial model. We summa- 
rize our principle results of the more massive LMC model 9 
as follows. 

(1) The formation of a global stellar bar cannot be com- 
pleted until T = —0.5 Gyr. This delayed formation of the 
stellar bar is due to the stronger self-gravity which prevents 
the Galactic tidal field from exciting the non-axisyrnmctric 
bar instability in the LMC disk. Only spiral arms are excited 
by the Galactic tidal field in the eaxly dynamical evolution of 
the LMC. As a natural result of the delayed bar formation, 
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Figure 24. The final distribution of old stars stripped from the 
LMC disk in an AitofT projection for the fiducial model (upper) 
and for the model 9 (lower). The preset locations (T = Gyr) of 
the LMC and the SMC are indicated by a large filled circle and 
a small filled one, respectively, whereas the locations of the LMC 
and the SMC ~ 9 Gyr ago (T = -8.8 Gyr) are indicated by a 
large open circle and a small open one, respectively. 



the bending/warping associated with the bar formation oc- 
curs around T = Gyr and is accordingly seen in the edge-on 
view of the LMC at T = Gyr (See Figure 22). 

(2) As shown in Figure 22, the disk cannot be thickened 
so much compared with the fiducial model owing to the rel- 
atively weak tidal perturbation from the Galaxy. Figure 23 
shows that the vertical velocity dispersion (ctz) consequently 
increases very slowly during the Galaxy-LMC-SMC interac- 
tion until the last collision between the LMC and the SMC 
at T = —0.2 Gyr. The final Uz is ~ 2.2 times larger than 
the initial value and thus is a factor of 1.4 smaller than 
that of the fiducial model with a smaller Mlmc (10^° M©). 
This result implies that kinematical evolution of the stellar 
disk of the LMC strongly depends on its mass during the 
Galaxy-LMC-SMC interaction. 

(3) The abrupt increases of CTz between —0.2 < T < 
Gyr seen in Figure 23 is due to the strong tidal perturbation 
caused by the LMC-SMC collision. There are no significant 
differences in the evolution of the half-mass radius (i?h) be- 
tween the isolated LMC model 10 and the model 9, because 
radial mass transfer, which drives a dramatic change of i?h, 
is not efficient both in the isolated model 10 and the model 
9. 

(4) The total mass of old stars stripped from the LMC 
disk is smaller in the more massive LMC model (model 9) 
than in the fiducial model, because the LMC is more strongly 
bounded by self-gravity in the more massive LMC model. As 
a result of this, the "tidal stream" composed of the stripped 
stars shown in Figure 24 is less clearly seen in the more 
massive model. This result implies that projected number 
density of the Galactic halo stars along the possible tidal 
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stream formed from stars stripped from the LMC stellar 
disk provides a valuable information of the LMC mass. 

(5) The formation efficiency of new field stars (esti- 
mated from the fraction of the mass of new stars to the 
initial gas mass) is similar between the more massive LMC 
model and the fiducial model (43 % and 47 %, respectively). 
All of the new GCs are formed only after T < — 3 Gyr and 
most of the newly formed GCs have the ages of < 2 Gyr in 
the model 9. 

3.2.2 Other parameters 

The parameter dependences of the present simulations are 
summarized as follows. 

(1) Final structural and morphological properties of the 
LMC disk and star formation histories of field stars and GCs 
do not depend strongly on 9 and (j) (i.e., inclination of the 
disk) for a reasonable set of these values. For example, mass 
fraction of gas converted into new field stars range from 40 
% to 47 % for the models 1, 3, 4, 5, 8 and all these models 
show a central stellar bar composed both of old field stars 
and new ones at T =0 Gyr. 

(2) The present results do not depend on the shape of 
the dark matter halo (models 5 and 6), essentially because 
the model parameters of the NFW halo are chosen such that 
Vm and the shape of the rotation curve arc similar to those 
of the fiducial model. The epoch of the commencement of 
the GC formation does not depend on the shape of the dark 
matter halo, which confirms that the most recent episode of 
globular cluster formation in the LMC is related to the com- 
mencement of strong tidal interactions between the LMC, 
the SMC and the Galaxy. 

(3) The star formation histories of field stars and GCs 
depend on the orbital evolution of the Clouds (model 11 
— 17), though the final barred morphology is not so differ- 
ent between the models. For example, the formation rate of 
field stars before T = —4 Gyr is significantly higher than 
that after T — —4 Gyr, for the model with the orbital type 
I (i.e., model 15), in which the LMC-SMC distance does not 
become so small compared with the fiducial model. As a re- 
sult of this, the age distribution of new field stars in this 
model 15 has no peak around T = —2 Gyr where the fidu- 
cial model shows its secondary peak in the age distribution: 
The age distribution of field stars provides a fossil record 
of the Galaxy-LMC-SMC interaction history that is totally 
determined by the orbits of these three. 

(4) The model without the SMC's tidal effect (model 
17) shows much less efficient formation of field star and 
GCs compared with the fiducial model with the SMC. For 
example, the mass fraction of gas converted into new field 
stars (new GCs) for the last ~ 9 Gyr in the model without 
the SMC is only 62 (17) % of the fiducial model with the 
SMC. This result strongly suggests that the tidal pertur- 
bation from the SMC plays an important role in the star 
formation history of the LMC, in particular, in the GC for- 
mation. 

(5) The dissipative dynamics of interstellar medium 
during the Galaxy-LMC-SMC interaction also controls the 
formation histories of field stars and GCs, though it can- 
not significantly change the age distributions of these stellar 
populations (models 1 and 18 — 21 with different /r). In 
particular, the formation efficiency of GCs (i.e., the total 



number of new GCs) strongly depends on /r in the sense 
that the efficiency is higher in the model with larger /r (i.e., 
more dissipative). 



4 DISCUSSION 

4.1 Born as a pair or difFerent entities ? 

It remains unclear whether the Clouds were born initially as 
a primordial pair of galaxies at the epoch of galaxy forma- 
tion (hereafter referred to as "the primordial binary galaxy 
model" just for convenience) or formed as different entities 
at different places and have only recently become a close pair 
for the first time ("the recent coupling model"). This ques- 
tion is the central core of any problems related to the long- 
term evolution of the Clouds (~ 10 Gyr), because strong 
dyuamical interaction between the Clouds significantly influ- 
ences not only star formation histories but also structure and 
kinematics of the Clouds if they axe strongly coupled. The 
primordial binary galax;y model has the following two advan- 
tages in explaining some physical properties of the Clouds. 
Firstly, previous numerical studies showed that structural 
and kinematical properties of the Magellanic stream are self- 
consistently explained by the models for which the Magel- 
lanic Clouds have remained in a binary state for the past 
15 Gyr (MF; GSF; GN; YN). Secondly, such a primordial 
binary galaxy model of the Clouds also naturally explains 
active star formation about 0.2 Gyr ago in the LMC disk 
(GSF) and structural properties of the stellar halo and the 
very recent star formation history of the SMC (GN; YN). 
MF discussed whether the binary Clouds could be formed 
from a single protogalaxy via tidal fission ~ 15 Gyr ago. 
Since these previous studies showed only binary galaxy mod- 
els that explain the above observational properties of the 
Cloud and the Magellanic Stream, it is not so clear whether 
other models for which the duration of the binary state of 
the Clouds does not exceed ~ 15 Gyr can equally explain 
physical properties. 

In calculating the orbital evolution of the Clouds, the 
above previous studies adopted the following two assump- 
tions that are not so consistent with the latest observations 
of LMC's structure and kinematics (e.g., vdMAHS): (1) the 
LMC has a mass of 2 x 10^° Mq and (2) dynamical fric- 
tion between the Clouds is negligible. The present study 
has demonstrated that if we adopt more reasonable assump- 
tions as to the above (1) and (2), it is not possible for the 
Clouds to remain in a binary state for the last 13 Gyr. It has 
therefore suggested that the LMC and the SMC could have 
dynamically coupled relatively recently (~ 4 Gyr ago) for 
the first time and thus were born not as a binary but as dif- 
ferent entities. These results never mean that the primordial 
binary galaxy models (e.g., MF) can be ruled out, because 
there are still some uncertainties in the mass estimation of 
the LMC by vdMAHS. It rather suggests that it depends 
on the model parameters. In particular, the LMC mass de- 
termines how long the Clouds can keep their binary state in 
their dynamical history: More precise mass estimation of the 
LMC is doubtlessly worthwhile for the better understanding 
of the orbital evolution of the Clouds. 

Formation histories of field stars and globular clusters 
and structural and kinematical properties in the LMC are 
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observed to be different from those in the SMC (e.g., van den 
Bergh 1981; 2000). For example, the "age gap" (i.e., only 
one cluster with the age ranging from ~ 13 Gyr to 3 Gyr) 
observed in LMC's globular clusters does not exist in SMC's 
globular clusters (e.g., Piatti ot al. 2002). Since strong tidal 
perturbation triggers the formation of globular clusters (e.g., 
Bekki et al. 2002), it is not clear why formation histories of 
star clusters are different between the two in the primordial 
binary model, in which the Clouds have been perturbing 
one another from their initial dynamical state until now. 
The observational fact that the age gap is seen only in the 
LMC can be understood in terms of the recent coupling 
model as follows. The LMC was formed as a relatively low 
surface brightness galaxy, being more distant (~ 150 kpc, 
corresponding to the apocenter of its early orbit) from the 
Galaxy so that the Galactic tidal field alone could not trigger 
cluster formation efficiently until it first encounters with the 
SMC. In contrast, the less massive SMC, which is therefore 
more susceptible to the Galactic tide, was born less distant 
(~ 100 kpc) from the Galaxy, and thus influenced by the 
Galaxy strongly enough to form globular clusters from the 
early evolutionary stage (several to lOGyrs ago). Thus the 
difference in cluster formation histories between the Clouds 
can be due to the differences in the birthplaces and initial 
masses between the two in the recent coupling model. 

Although Magellanic-type galaxies that appear to be 
pairs of galaxies are not rare (e.g.. Freeman 1984), it is not 
so obscrvationally clear whether these apparently pairs of 
galaxies were formed initially as pairs at the epoch of galaxy 
formation or have only recently become pairs owing to recent 
tidal capture or interaction of galaxies (e.g., Helou 1984). It 
is accordingly difficult for us to derive some hints on the 
above question as to the binary state of the Clouds from 
observations on other Magellanic- type pair. Currently avail- 
able proper motion data for the Clouds has the accuracy of 
~ 3 mas yr^^ (corresponding to ~ 700 km s~^ error of a 
single star for the distance to the Clouds) so that we cannot 
determine the duration of the LMC-SMC binary state di- 
rectly from observational data on the proper motion and the 
radial velocities of the Clouds (e.g., Westerlund 1997). Fu- 
ture astrometric missions with the ~ 10 /las accuracy (cor- 
responding to a few km s~^) will allow us to derive an 
unambiguous answer for the problem as to the binary state 
of the Clouds. 

4.2 Stellar halo formation in the LMC and 
Magellanic-type dwarfs 

Old, metal-poor stellar halo populations of a galaxy have 
long been considered to be "fossil records" which contain 
valuable information on dynamical and chemical evolution 
of the galaxy (e.g., Eggen, Lynden-Bell, & Sandage 1962). 
Recent numerical simulations based on the currently favored 
cold dark matter (CDM) theory of galaxy formation have 
demonstrated that basic physical processes involved in the 
formation of the stellar halo are described by both dissipa- 
tive and dissipationless merging of subgalactic clumps and 
their resultant tidal disruption in the course of gravitational 
contraction of the Galaxy at high rcdshift (Bekki & Chiba 
2000, 2001). Previous studies also suggested that (1) merg- 
ing of subgalactic clumps is an essentially important process 
for the stellar halo formation in disk galaxies like the Galaxy 



and (2) the physical properties of the halos depend on the de- 
tails of the merging processes of subgalactic clumps (Bekki 
& Chiba 2000, 2001). It is, however, unclear (1) whether 
stellar halos exist in less luminous disk (or irregular/dwarf) 
galaxies like the LMC and (2) how they are formed in the 
course of their formation. 

The present study has demonstrated that (1) the stel- 
lar halo is formed from redistribution (in space) of stars 
initially within the outer part of the disk in the LMC even 
if the LMC has initially no stellar halo, (2) the developed 
stellar halo contains some fraction of younger stars, and (3) 
the spatial distribution of the outer stellar halo is not ho- 
mogeneous. Accordingly the present study suggests that the 
formation process of the LMC's stellar halo differs from the 
Galactic one in that it is not associated with any merging 
of subgalactic clumps. The above result (1) furthermore im- 
plies that (i) the LMC could either have no or little amount 
of stellar mass in its halo at the epoch of its formation, even 
if it is now observed to have the kinematically hot stellar 
halo (Minniti et al. 2003; Alves 2004), and (ii) it acquired 
the substantial mass of the halo owing to the strong tidal 
interaction between the Galaxy and the SMC (This could 
be true for the SMC, which could suffer more severe tidal 
perturbation from the Galaxy and the LMC than the LMC 
could). The present results thus raise the following two is- 
sues as to the stellar halo formation in less luminous disk 
galaxies: (1) whether less luminous disk galaxies are formed 
initially without old stellar halos and (2) if so, why no stel- 
lar halos are formed in such galEixies at the epoch of galaxy 
formation. 

Recently the physical properties of old stellar halos in 
such less luminous galaxies have been investigated, in par- 
ticular, in dwarf galaxies in the Local Group (e.g., Demers et 
al. 2003). One of the intriguing results is that NGC 3109, a 
Magellanic-type dwarf on the outskirts of the Local Group, 
contains carbon stars nearly exclusively in and near its disk 
component, whereas NGC 6822, a galaxy with the same mor- 
phological type, has an extended intermediate-age halo as 
well as an old halo. This suggests that the formation of an 
old stellar halo is strongly suppressed in some of Magellanic- 
type dwarfs like NGC 3109. The following two mechanisms 
are considered for suppressing the formation of a stellar halo. 
One is that while the formation of disk galaxies like NGC 
3109 may accompany merging/accretion incidents of sub- 
galactic clumps, if such clumps were totally gaseous with- 
out containing old stars, the accretion/merging of gaseous 
clumps would not leave a diffuse halo component, which 
originates from tidal disruption of pre-existing stars. The 
other is that the initial gas distribution in these galaxies 
were so diffuse that star formation could happen only af- 
ter the settlement of gas onto the disk plane, whereby the 
regions with high gas density emerged. 

If we adopt the currently favored theory based on hier- 
archical assembly of CDM (White & Rees 1978), the above 
second scenario seems less likely, because the CDM model 
predicts higher overdensities for galaxies embedded in less 
massive dark halos. Then, if the first scenario is the case, 
what is the most likely mechanism for the suppression of star 
formation within subgalactic clumps which end up with less 
luminous (Magellanic-type dwarf) disk galaxies? One pos- 
sible mechanism is that thermal and/or kinematic feedback 
supplied by supernovae significantly suppresses star forma- 



© 1994 RAS, MNRAS 000, 1-?? 



Formation of the LMC 21 



tion in subgalactic clumps which end up with less luminous 
disk galeixies (e.g., Dekel & Silk 1986). In this mechanism, 
less massive clumps are susceptible to supernova explosions 
owing to their shallower gravitational potential. Alterna- 
tively, the UV background in the Universe suppresses the 
formation of dwarf galaxies via photoionization effect (e.g.. 
Bullock, Kravtsov & Weinberg 2000), in such a way that 
less massive galaxies with lower virial temperature may be 
more affected. We note that these mechanisms are also clues 
to solving the problem of overabundance in the number of 
CDM subhalos in the Local Group (Klypin et al. 1999). Bet- 
ter understanding of such suppression effects of galaxy for- 
mation may also resolve the current issue of stellar halos in 
less luminous disk galaxies. 

The present study predicts that if the major compo- 
nent of the stellar halo in a less luminous (Magellanic-type 
dwarf) disk galaxy is formed by tidal interaction with other 
more luminous galaxies, the disk galaxy has (1) relatively 
young stellar halo populations, most of which come from 
the outer disk stars formed before the interaction and (2) the 
disk has a thick disk formed thorough the disk heating by 
tidal interaction. Therefore, future observations on (1) age 
and metallicity distributions of halo stars in Magcllanic-typc 
dwarf disk galaxies and (2) statistical correlations between 
the presence of the outer stellar halos and that of the faint 
thick disk components will help us to determine whether 
the stellar halos in such dwarfs are formed from tidal galaxy 
interaction rather than from primordial merging/accretion 
of subgalactic clumps. Heidmann et al. (1972) revealed that 
the intrinsic flattening in disk galaxies decreases (i.e., less 
flattened) abruptly from Sm to Im Hubble types. If this less 
flattened nature is due to disk heating in these dwarf irreg- 
ulars, it is an observationally interesting question (related 
to the stellar halo formation via tidal interaction) whether 
spherical stellar halos are more likely to be observed in Im 
rather than Sd galaxies. In the stellar halo formation sce- 
nario via tidal interaction, the age distribution of halo stars 
in a galaxy depends strongly on when the galaxy interacted 
with other galaxy. Therefore, future observations on the age 
distribution of halo stars will also provide valuable informa- 
tion on the past interaction history of the galaxy. 

4.3 Origin of the LMC's stellar bar 

Recently several observational studies have attempted to de- 
rive age and metallicity distributions of stellar populations 
in the bar region of the LMC in order to constrain the star 
formation history (SFH) in the bar (Elson et al. 1997; Ardel- 
berg et al. 1997; Holzman et al. 1999; Olsen 1999; Smecker- 
Hane et al. 2002). There however exists some discrepancy 
in the results of the SFH of the bar between different obser- 
vations, possibly because authors investigated SFHs of dif- 
ferent regions within the bar using different number of stars 
analyzed (e.g., Smecker-Hane et al. 2002). For example, El- 
son et al. (1997) investigated photometric properties of ~ 
15800 stars obtained by the HST for the inner LMC disk 
and found a possible evidence of a later starburst around 
1 Gyr ago which may be responsible for the bar formation 
in the LMC disk. Smecker-Hane et al. (2002) revealed that 
star formation of the dominant populations in the LMC bar 
occurred from 4 to 6 and 1 to 2 Gyr ago. Olsen (1999) sug- 
gested that the LMC bar region appears to have high levels 



of star formation activity as long as 5 — 8 Gyr ago: The LMC 
bar is dominated by old stellar populations. 

The present study has demonstrated that (1) a large 
fraction of stars (up to 50 %) in the bar are formed dur- 
ing the strong tidal interaction between the Clouds and the 
Galaxy and have relatively younger ages and (2) there can 
exist a steep age gradient of stellar populations along with 
and perpendicular to the bar in the sense that the outer 
regions of the bar contain only a smaller fraction of young 
stars. We therefore suggest that the formation of the ob- 
served young stellar populations in the LMC bar region (e.g., 
Elson et al. 1997; Smecker-Hane et al. 2002) is closely as- 
sociated with the efficient star formation within the bar for 
the last several Gyr (in particular, ~ 2 Gyr ago). We also 
suggest that the observed discrepancy in the SFH of the bar 
could be due partly to the radial gradient of ages of stellar 
populations within the bar: The meiss fraction of young stars 
and the mean age of stellar populations for a target field of 
the LMC bar in previous observations depend strongly on 
the distance of the field from the center of the LMC. 

The present results imply that the LMC bar was formed 
not spontaneously from global bar instability in the early 
evolution stage of the LMC but from tidal perturbation by 
the Galaxy and the SMC. Several numerical studies have 
already shown that stellar bar can be formed via tidal in- 
teraction of disk galaxies for variously different parameters 
of galaxy interaction (e.g., Noguchi 1987; Byrd & Valtonen 
1990). The mass fraction of stars within a disk embedded 
by a massive dark matter halo must be at least larger than 
0.4 — 0.5 so that the bar is spontaneously formed from global 
bar instability (e.g., Sellwood & Caxlberg 1984). The total 
visible mass of stars of the LMC with Ly = 3.0 x 10^ Lq is 
~ 2.7 X lO'' Mq for M/Lv = 0.9 ± 0.2) whereas the total 
dynamical mass of the LMC is (8.7 ± 4.3) x 10^ M© within 
8.9 kpc (vdMAHS). Therefore, the spontaneous bar forma- 
tion in the LMC disk is not likely to occur and thus the bar is 
likely to have formed relatively recently from external tidal 
perturbation. Using numerical simulations, Noguchi (1996) 
demonstrated that stellar bars formed from external tidal 
perturbation ( "tidal bars" ) have a relatively flat density pro- 
file along the major axis of the bars with "shoulders" (abrupt 
steepening of the gradient) at the bar ends. We thus suggest 
that future observational studies on the radial density pro- 
file for the young stellar populations of the LMC confirm 
the bar's radial profiles characteristic of the tidal bars, if 
the LMC bar was formed from tidal interaction with the 
Galaxy and the SMC relatively recently. 

4.4 The Age gap problem 

Precise estimation of an age of each individual star cluster 
in the LMC leads to the determination of the cluster age dis- 
tribution and thus to the better understanding of the star 
formation history of the LMC (e.g., Searle et al. 1980; Hodge 
1983, 1988; Mateo 1988). Differences in spatial distributions 
between clusters with different ages provide some informa- 
tion on the spatial variation of the star formation history of 
the LMC (e.g., van den Bergh 1981). The age distribution 
of the LMC clusters shows a gap extending from 13 to 3 Gyr 
with only one cluster (ESO 121-SC03) within this gap, which 
is not seen in the SMC clusters (Jensen et al. 1988; Da Costa 
1991; Geisler et al. 1997; Rich et al. 2001; Piatti et al. 2002). 
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The following three possible scenarios are proposed for ex- 
plaining the above "age gap" in the LMC clusters. First is 
that star cluster formation after the initial formation of old 
globular clusters ~ 13 Gyr ago (at the epoch of the LMC 
formation) had been suspended until very recently ^ 3 Gyr 
ago. The second scenario is that although cluster formation 
has ceaselessly continued until now, only star clusters with 
ages ranging from 13 to 3 Gyr are tidally stripped. Third 
is that star clusters with ages between 3 and 13 Gyr were 
preferentially destroyed by the LMC tidal field to become 
field stars. 

For the above second scenario to be viable, the star 
clusters with ages between 3 and 13 Gyr should be formed 
preferentially in the outer LMC's halo, where tidal stripping 
of the clusters by the Galaxy is very efficient in this scenario. 
Given the possible observational evidence that both young 
(< 3 Gyr old) and old (~ 13 Gyr) clusters show disky distri- 
butions (e.g., Schommcr ct al. 1991; van den Bcrgh 2000), 
it is unclear why only clusters with ages between 3 and 13 
Gyr are formed in the outer halo region of the LMC. There- 
fore the second scenario is regarded as rather less likely one. 
Regarding the third scenario, it could be possible for clus- 
ters with ages between 3 and 13 Gyr to be preferentially 
destroyed, only if they have typically lower densities and 
masses compared with other LMC clusters. Since no previ- 
ous theoretical models predicted age dependences of struc- 
tural properties of globular clusters (e.g., Harris 1991 and 
reference therein), the third scenario is equally less likely 
one. 

Thus, if the first scenario is only a reasonable one, the 
essence of the problem related to the above "ago gap" in 
the LMC clusters is what mechanism is responsible for the 
abrupt reactivation of cluster formation in the LMC ^ 3 
Gyr. Using numerical simulations, Bekki et al. (2004b) first 
discussed this problem in the context of mutual tidal in- 
teraction between the Clouds and the Galaxy. Bekki et al. 
(2004b) and the present study have proposed that the epoch 
of reactivation of cluster formation corresponds to the com- 
mencement of strong tidal interaction between the LMC and 
the SMC, which disturbs the LMC gas disk, enhances cloud- 
cloud collision rate, and consequently triggers cluster for- 
mation. In this scenario, the tidal interaction between the 
LMC and the Galaxy alone cannot increase so dramatically 
the number of cloud-cloud collisions leading to cluster for- 
mation between 13 and 3 Gyr ago. 

The above scenario could be just one of promising sce- 
narios explaining the origin of the age gap, we accordingly 
point out two possible alternative scenarios below. Byrd et 
al. (1994) numerically investigated the orbital evolution of 
the Clouds and Leo I and thereby proposed a scenario that 
the Cloud left M31 ~ 10 Gyr ago and were tidally cap- 
tured by the Galaxy several Gyr ago. Although they did 
not discuss their results in terms of the age gap of the LMC 
clusters, it is not unreasonable to expect that the cluster for- 
mation could be suddenly triggered by the strong Galactic 
tidal force when the LMC first experienced the pericenter 
passage with respect to the Galaxy. One of the alternative 
scenarios is thus that the origin of the age gap is closely asso- 
ciated with the first pericenter passage of the LMC that once 
belonged to M31. In this scenario, Byrd et al. (1994) showed 
that the epoch of the first pericenter passage is ranging from 
4.6 Gyr to 12.2 Gyr ago for a relatively narrow parameter 



space of the orbital evolution. Therefore it remains loss clear 
whether LMC's first pericenter passage is not ~ 6 Gyr ago 
but ~ 3 Gyr ago. Numerical studies with more variously dif- 
ferent yet reasonable initial orbital parameters of the LMC 
and with more realistic mass models of the Local Group will 
thus assess the viability of this scenario. 

The other alternative scenario is that ~ 3 Gyr ago cor- 
responds to the epoch when the LMC's disk gas begins to 
interact with the Galactic halo plasma and form clusters 
owing to bow-shocked induced star formation in the LMC 
disk, de Boer et al. (1998) pointed out that a shock induced 
by the ram pressure of the halo plasma can induce star for- 
mation in the LMC for a reasonable set of parameters of 
halo gaseous density and temperature and the relative ve- 
locity of the LMC with respect to the Galaxy. Hydrody- 
namical simulations demonstrated that ram pressure of the 
intracluster/intragroup medium strongly compresses a self- 
gravitating gas cloud within a short time scale (~ lO'^ yr), 
dramatically increasing the central gas density, and conse- 
quently causing efficient formation of a compact star cluster 
within the cloud (Bekki & Couch 2003) . Although these two 
works suggest that cluster formation via ram pressure of the 
Galactic halo plasma is possible, it is not clear why cluster 
formation via ram pressure becomes possible for the first 
time only ^ 3 Gyr ago, given the fact that the LMC had 
experienced several pericenter passages before 3 Gyr ago. 

One of the possible reason for the sudden interaction 
between the LMC gas and the Galactic halo gas ~ 3 Gyr 
ago is that dynamical friction from the Galactic dark mat- 
ter halo causes the orbital decay of the LMC so slowly that 
the LMC cannot approach so closely the outer edge of the 
Galactic hot plasma until recently (~ 3 Gyr ago). There 
could be more reasonable scenarios for the age gap problem 
other than the three discussed above. The age gap reflects 
the complicated interaction history between the Clouds and 
the Galaxy, which depends almost exclusively the details of 
the orbital evolution of the Clouds. This emphasizes the im- 
portance of proper motion measurement with the accuracy 
of an order of pas for the Clouds in clarifying the origin of 
the age gap of the LMC clusters. 



5 CONCLUSIONS 

We have performed numerical simulations for the dynami- 
cal and chemical evolution of the LMC interacting with the 
Galaxy and the SMC for the last 9 Gyr. The main results 
are summarized as follows. 

(1) Tidal interaction between the Clouds and the 
Galaxy plays a major role not only in the morphological 
transformation of the LMC disk but also in the formation 
history of field stars and star clusters. The interaction trans- 
forms the initially thin, non-barred LMC disk into the three 
different components; the thick disk, bar, and kinematically 
hot stellar halo. The central bar formed during the tidal 
interaction is composed both of old field stars and newly 
formed ones with each fraction being equal in the innermost 
part. The final thick disk has the velocity dispersion of ~ 30 
km s^ and shows rotationally supported kinematics with 
Vm/ao ^ 2.3. The outer stellar disk (i?L > 5 kpc) surround- 
ing the central bar is highly elliptic resulting from the tidal 
interaction between the Galaxy and the SMC. 
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(2) The stellar halo is formed during the interaction as 
a result of redistribution of stars initially within the outer 
part of the thin LMC disk. The stellar halo is thus com- 
posed mainly of old stars originating from the outer part of 
the initially thin LMC disk. The outer halo shows velocity 
dispersion of ~ 40 km s~^ at the distance of 7.5 kpc from 
the LMC center and has somewhat inhomogeneous distri- 
bution of stars. The stellar halo contains relatively young, 
metal-rich stars with the mass fraction of 2 % in the halo. 
Therefore the MDF of the halo is determined by the old 
stars and shows a peak around [Fe/H] ~ —1.3 for a reason- 
able set of parameters of the initial MDF in the stellar disk 
of the LMC. 

(3) Star formation rate in the LMC disk is moderately 
and repeatedly enhanced owing to the repetitive interac- 
tion between the Clouds and the Galaxy. The star forma- 
tion rate increases from ~ 0.1 Mq yr~^ to ~ 0.4 Mq yr~^ 
at the first periccnter passage of the LMC with respect to 
the Galaxy about 7 Gyr ago. The star formation rate also 
becomes moderately high when the Cloud begins to inter- 
act violently (with the LMC-SMC pericenter of less than 10 
kpc) about ~ 3.5 and 2 Gyr ago. Most of the new stars (~ 
90 %) arc formed within the central 3 kpc, in particular, 
within the bar for the Icist 9 Gyr. Consequently, the half 
mass radius is different by a factor of 2.3 between old field 
stars and newly formed ones. These structural differences 
between field stars with difference ages are characteristic of 
the LMC disk under tidal interaction with the Galaxy and 
the SMC. 

(4) Efficient GC formation does not occur until the 

LMC starts interacting violently and closely with the SMC 
(~ 4 Gyrs ago). This is due to the fact that cloud-cloud 
collisions with moderately high relative speed (30 < Kci < 
100 km s~^) and with small impact parameter (b < 0.25) 
required for GC formation occurs the most frequently when 
both the SMC and the Galaxy dynamically influence the 
LMC's disk strongly. The newly formed GC system has a 
disky distribution with rotational kinematics and its mean 
metallicity is ~ 1.2 higher than that of new field stars be- 
cause of the pre-enrichment by the formation of field stars 
prior to cluster formation. 

(5) About 15(20) % of the field stars (gas) initially 
within the LMC disk are tidally stripped to form a great 
stellar (gaseous) circle of a relic stream around the Galaxy 
during the last 9 Gyr evolution of the LMC. The great stellar 
circle shows inhomogeneity in some parts and is composed 
only of metal-poor old stars. The unique distributions of dis- 
tance and radial velocity in the tidal stream may well enable 
xis to pick out the stream among the Galactic halo stars. The 
stellar total mass of the tidal stream depends on the initial 
mass of the LMC (i.e., smaller for the larger LMC mass), so 
that the stellar number density along the stream provides 
valuable information on the LMC mass. 

(6) The LMC evolution depends on its initial mass and 
orbit with respect to the Galaxy and the SMC. In particular, 
the epoch of the bar and the thick disk formation is deter- 
mined by the LMC mass in such a way that the stellar bar 
and the thick disk are formed later in the model with a larger 
LMC mass. The mass fraction of the stellar halo is smaller 
for the model with a larger LMC mass. These are essen- 
tially because the LMC with a larger mass is more strongly 
bounded by its self-gravity so that the tidal perturbation 



from the Galaxy and the SMC does not so significantly in- 
fluence the dynamical evolution of the LMC. 

Based on these results, we have discussed the origin of 
the stellar halo in less luminous late-type galaxies such as 
the Clouds, the formation of LMC's stellar bar and thick 
disk, the origin of the age gap of the LMC globular clus- 
ters, and the difference in formation histories of field stars 
and globular clusters between the LMC and the SMC. It has 
been also pointed out that future proper motion measure- 
ments of the Clouds with ~ 10 /lac accuracy (Perryman et 
al. 2001) to estimate their past 3D orbits in an unprecedent- 
edly precise manner will provide us invaluable information 
on the complicated interplay between dynamical evolution 
of the LMC and its star formation history. 
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APPENDIX A: STATISTICS ON THE 
DURATION OF THE LMC-SMC BINARY 
STATUS 

As suggested by Bokki et al. (2004b) and the present study, 
the epoch when the LMC and the SMC become dynamically 
coupled (i.e., when the pericenter distance of the LMC-SMC 
orbit is as small as 10 kpc) is a critical moment for the 
LMC, because the combined tidal effect of the Galaxy and 
the SMC starts influencing significantly the LMC evolution 
after the dynamical coupling. It is thus important to in- 
vestigate when they are the most likely to become coupled 
during their dynamical evolution. Although the proper mo- 
tion of the Clouds have been already derived by several au- 
thors (e.g., KB; vdMAHS), the measurement error in these 
observational studies is so large that precise estimation of 
current velocity components {Ul,Vl,,Wl,) and {Us,Vs,Ws) 
within an error of ~ a few km s~^ cannot be made directly 
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Table Al. Orbit models 



orbit model 
A 
B 
C 
D 



-^bin,l(%) 




16.7 
43.8 



fbin.2(%) 

5.5 X 10-4 
1.1 X 10^2 
1.1 X 10-'^ 

4.6 X 10-2 



tbin,i(Gyr) 
4.29 
4.42 
6.25 



*bin,2(Gyr) 
1.15 
1.10 
1.23 
1.26 



0.8 
0.6 
0.4 
0.2 
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Figure Al. The distributions of t^in duration of binary status of 
the Clouds for the models with Mlmc = lO"*^" Mq (upper) and 
those with Mlmc = 2 x lO'^'^ Mq (lower), t^in represents the 
duration of binary status of the Clouds and the definition of the 
binary status is given in the main text. For comparison, the results 
are shown for models with (solid) and without (dotted) dynam- 
ical friction between the Clouds in each panel. For convenience, 
the normalized number of models is given for each bin of the 
LMC-SMC binary duration. Here only the models in which each 
component of the current velocities of the Clouds (((/l, Vl, Wl) 
and ([^Si ^S) M^s)) within ± 10 km of the corresponding 
component in the model. 



from the observations. As shown in previous studies on the 
orbital evolution of the Clouds (e.g., MF), only a velocity 
difference of ~ 5 km s~^ can cause a significant difference 
in the orbital evolution and thus in the binary status of the 
Clouds. Thus currently available data alone do not allow us 
to make a robust conclusion on the duration of the LMC- 
SMC binary status. 

Therefore we here make a statistical argument on the 
duration of the LMC-SMC binary status by investigating 
every possible orbits of the Clouds for a set of reasonable 
parameters of gravitational potential for the Galaxy and the 
Clouds. Using the backward integration scheme (MF) and 
the models for gravitational potential and dynamical friction 
adopted in this study for the Galaxy and the Clouds (i.e., the 
equation (1) - (5)), wc first calculate the orbital evolution of 
the Clouds for each model with a given set of parameters for 
their initial velocities. Then we estimate the epoch when the 
LMC-SMC distance is for the first time larger than 50 kpc 
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Figure A2. The same as Figure Al but for the models with the 
current velocities of the Clouds similar to those obtained from 
Hipparcos data (KB). 



that is adopted by previous numerical studies (MF; GSF) 
for the threshold of the binary status. 

Initial velocity components of Uh and Vl (f^s and Vs) 
for the LMC (the SMC) range from -320 km s"^ to 320 km 
s~^ in the orbital investigation. Wl (Ws) is derived from 
given Uj, and Vl {Us and Vs) in the LMC (the SMC) so 
that the radial velocity of the LMC (the SMC) is consistent 
with the observed radial velocity of the LMC (the SMC). Wo 
survey the LMC-SMC binary orbit every 3.2 km s^^ in the 
range of {Uh,VL,Wh) -space {{Us,Vs,Ws)-space) and thus 
the total number of the orbits investigated in this study is 
an order of (10^ x 10^) x (10^ x 10^) = 10^ Such a numerous 
number of orbital investigation enable us to provide a sta- 
tistical argument on the duration of the LMC-SMC binary 
status. 

We adopt the four different orbital models with differ- 
ent Mlmc and with or without dynamical friction: Model 

01 with Mlmc = 10^° Mq and with dynamical friction, 

02 with Mlmc = 10^° Mq and without dynamical fric- 
tion, 03 with Mlmc = 2 x 10^° Mq and with dynamical 
friction, and 04 with Mlmc = 2 x 10^° Mq and without 
dynamical friction. We investigate not only the models with 
Mlmc consistent with observations (i.e., Ol and 02) but 
also those with Mlmc significantly larger than the observed 
one (vdMAHS), because we intend to compare the present 
results with previous ones (e.g., GN) in which Mlmc is as- 
sumed to be 2 X 10^° Mq. 

Based on the above models, we search for the models 
which meat the following two requirements/conditions: (CI) 
The current velocities are broadly consistent with those de- 
rived from the latest observations, and (C2) the Clouds can 
keep its binary status for the Hubble time (~ 13 Gyr). We 
first determine the models satisfying the above condition CI 
by selecting the models with each of the current velocities 
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(e.g., Ul) within an error of 10 km s^^ of the observationally 
inferred one among all models. We then search for the mod- 
els satisfying the above condition C2 among those satisfying 
the above condition CI based on the backward integration 
of the Clouds' orbits. In the first selection process, wc check 
whether the current velocities in each model are consistent 
either with those adopted by GN or with those by KB. 

We mainly investigate (1) number fraction of orbits 
(Jbin) satisfying the above two conditions among all pos- 
sible ones and (2) mean duration of the LMC-SMC binary 
status (thin) for each of the four models, Ol, 02, 03, and 04. 
The number fraction of orbital models in which the above 
condition CI is satisfied and the current velocities arc within 
an error of 10 km s^^ of those by GN (KB) is represented by 
^bin,i (^bin,2)- The mean duration of the LMC-SMC binary 
status for orbits with the current velocities within an error 
of 10 km s~^ of those by GN (KB) is represented by tbin.i 
(ibin,2)- The Table Al summarizes the results of the orbital 
investigation; Model number (column 1), Fbin.i (2), Fbin,2 
(3), tbin.i (4), and tbin,2 (5). 

Figure Al shows the distribution of tbin.i of orbits in 
which the current velocities are within an error of 10 km 
s~^ of those by GN for the 01 — 04 models. It is clear from 
this figure that (1) the Clouds can keep their binary state for 
only less than 7 Gyr (i.e., the Clouds becomes disintegrated 
until T = -7 Gyr) in the models 01 and 02, (2) they are 
the most likely to keep their binary status for ~ 4 Gyr in 
the models 01 and 02, (3) tbin,i is the most likely to be ~ 2 
Gyr for the 03 model, and (4) the Clouds are the most likely 
to be able to keep their binary status for the Hubble time 
only if the LMC has a large mass of 2. O x M© without 
dynamical friction between the Clouds. These results imply 
that the Clouds are very hard to keep their binary status for 
the Hubble time if a reasonable sets of assumptions (Mlmc 
= 10^" Mq and inclusion of the dynamical friction between 
the Clouds) are made for the orbital calculations. 

Figure A2 shows the distribution of tbin,2 for orbital 
models in which the current velocities are within an error of 
10 km s~^ of those by KB for the Ol — 04 models. As shown 
in this Figure A2, the Clouds are the most likely to be able 
to keep their binary status for less than 1 Gyr, irrespectively 
of the LMC mass and whether or not the dynamical friction 
between the Clouds is included in the orbital calculations. 
The more massive LMC has a higher probability of keep- 
ing the LMC-SMC binary status (i.e., larger Fhin,2) for the 
Hubble time in these models with and without dynamical 
friction (See the Table Al). The results shown in Figure Al 
and A2 thus strongly suggest that the probability of Clouds 
keeping their binary status is very low. 
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